Chapter 1 



Introduction 



1.1 What is a Plasma? 
1.1.1 An ionized gas 

A plasma is a gas in which an important fraction of the atoms is ionized, so that the electrons 
and ions are separately free. 

When does this ionization occur? When the temperature is hot enough. 
Balance between collisional ionization and recombination: 



Ionization 




Figure 1.1: Ionization and Recombination 

Ionization has a threshold energy. Recombination has not but is much less probable. 
Threshold is ionization energy (13.6eV, H). \% 

Integral over Maxwellian distribution gives rate coefficients (reaction rates). Because of 
the tail of the Maxwellian distribution, the ionization rate extends below T — Xi- And in 
equilibrium, when 

^ions _ < CTjV > 
^-neutrals ^ 0~ r V i> 
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Figure 1.2: Ionization and radiative recombination rate coefficients for atomic hydrogen 

the percentage of ions is large (~ 100%) if electron temperature: T e > Xi/10. e.g. Hydrogen 
is ionized for T e > leV (ll,600°k). At room temp r ionization is negligible. 

For dissociation and ionization balance figure see e.g. Delcroix Plasma Physics Wiley (1965) 
figure 1A.5, page 25. 



1.1.2 Plasmas are Quasi-Neutral 

If a gas of electrons and ions (singly charged) has unequal numbers, there will be a net charge 
density, p. 

p = n e (-e) + rii(+e) = e(rii - n e ) (1.2) 
This will give rise to an electric field via 



V.E = ^ = -{m-n e ) 



(1.3) 



Example: Slab. 



dE 

dx 

-> E 



x 



= P- 



(1.4) 
(1.5) 



7 



+ 

/ 1 

+/ + 
+ 



Figure 1.3: Charged slab 

This results in a force on the charges tending to expel whichever species is in excess. That 
is, if rii > n e , the E field causes n« to decrease, n e to increase tending to reduce the charge. 

This restoring force is enormous] 
Example 

Consider T e = leV ', n e = 10 19 m~ 3 (a modest plasma; c.f. density of atmosphere n mo i CC uies ~ 
3 x 10 25 m~ 3 ). Suppose there is a small difference in ion and electron densities An = (n { — n e ) 



so p = An e 
Then the force per unit volume at distance x is 

F e = pE = p 2 - = {An e) 2 — 

Take An/n e = 1% ,x = 0.10m. 

F e = (10 17 x 1.6 x 10" 19 ) 2 0.1/8.8 x 10" 12 = 3 x 10 6 A^.m" 3 
Compare with this the pressure force per unit volume ~ p/x : p ~ n e T e (+niTi) 

F p ~ 10 19 x 1.6 x 10" 19 /0.1 = 16iVm- 3 



1.6) 



(1-7) 



(1.8) 



(1.9) 



Electrostatic force >> Kinetic Pressure Force. 

This is one aspect of the fact that, because of being ionized, plasmas exhibit all sorts of col- 
lective behavior, different from neutral gases, mediated by the long distance electromagnetic 
forces E, B. 

Another example (related) is that of longitudinal waves. In a normal gas, sound waves are 
propagated via the intermolecular action of collisions. In a plasma, waves can propagate 
when collisions are negligible because of the coulomb interaction of the particles. 
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1.2 Plasma Shielding 

1.2.1 Elementary Derivation of the Boltzmann Distribution 

Basic principle of Statistical Mechanics: 

Thermal Equilibrium <-> Most Probable State i.e. State with large number of possible ar- 
rangements of micro-states. 



Figure 1.4: Statistical Systems in Thermal Contact 

Consider two weakly coupled systems Si, S 2 with energies Ei,E 2 . Let gi, g 2 be the number of 
microscopic states which give rise to these energies, for each system. Then the total number 
of micro-states of the combined system is (assuming states are independent) 

9 = 9x92 (1-10) 

If the total energy of combined system is fixed E x + E 2 = E t then this can be written as a 
function of E\. 

g = g 1 (E 1 )g 2 (E t -E 1 ) (1.11) 
, dg dg x dg 2 
and dEl = dE 92 ~ 9l dE ■ (L12) 

The most probable state is that for which jjjr = i.e. 

1 dgx 1 dg 2 d d 

= or — moi = — mo? (l-lo) 

gi dE g 2 dE dE yi dE y2 K J 

Thus, in equilibrium, states in thermal contact have equal values of ^lng. 
One defines a = In g as the Entropy. 
And [^lng] -1 = T the Temperature. 

Now suppose that we want to know the relative probability of 2 mzcro-states of system 1 in 
equilibrium. There are, in all, g± of these states, for each specific E\ but we want to know 
how many states of the combined system correspond to a single microstate of Si. 

Obviously that is just equal to the number of states of system 2. So, denoting the two values 
of the energies of Si for the two microstates we are comparing by Ea, Eb the ratio of the 
number of combined system states for Sia and Sib is 

f# ~ t A \ = exp[(7(^ - E A ) - a(E t - E B )\ (1.14) 
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Now we suppose that system S 2 is large compared with Si so that Ea and E B represent very 
small changes in SVs energy, and we can Taylor expand 



g 2 (E t - E A ) 

— — — — ~ exp 

g 2 {E t - Ea) 



da da 
~ EA dE +EB dE 



;i.i5) 



Thus we have shown that the ratio of the probability of a system (S\) being in any two 
micro-states A, B is simply 

" -{Ea — Eb) 
T 



exp 



;i.i6) 



when in equilibrium with a (large) thermal "reservoir" . This is the well-known "Boltzmann 
factor" . 

You may notice that Boltzmann's constant is absent from this formula. That is because of 
using natural thermodynamic units for entropy (dimensionless) and temperature (energy). 

Boltzmann's constant is simply a conversion factor between the natural units of temperature 
(energy, e.g. Joules) and (e.g.) degrees Kelvin. Kelvins are based on °C which arbitrarily 
choose melting and boiling points of water and divide into 100. 

Plasma physics is done almost always using energy units for temperature. Because Joules 
are very large, usually electron-volts (eV) are used. 

leV = 11600K = 1.6 x 10" 19 Joules. (1.17) 

One consequence of our Botzmann factor is that a gas of moving particles whose energy is 
7,mv 2 adopts the Maxwell-Boltzmann (Maxwellian) distribution of velocities cx exp[— ^-]. 

1.2.2 Plasma Density in Electrostatic Potential 

When there is a varying potential, 0, the densities of electrons (and ions) is affected by it. 
If electrons are in thermal equilibrium, they will adopt a Boltzmann distribution of density 

,ecj) 

n e oc exp(— ) . (1.18) 
This is because each electron, regardless of velocity possesses a potential energy — e<p. 



Consequence is that (fig 1.5) a self-consistent loop of dependencies occurs. 

This is one elementary example of the general principle of plasmas requiring a self-consistent 
solution of Maxwell's equations of electrodynamics plus the particle dynamics of the plasma. 

1.2.3 Debye Shielding 

A slightly different approach to discussing quasi-neutrality leads to the important quantity 
called the Debye Length. 

Suppose we put a plane grid into a plasma, held at a certain potential, (j) g . 
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Figure 1.5: Self-consistent loop of dependencies 




Figure 1.6: Shielding of fields from a 1-D grid. 



Then, unlike the vacuum case, the perturbation to the potential falls off rather rapidly into 
the plasma. We can show this as follows. The important equations are: 

d 2 (f) e 

Poisson's Equation V 2 = -—^ = (n* — n e ) (1-19) 

dx 2 e 

Electron Density n e = exp(e0/T e ). (1-20) 

[This is a Boltzmann factor; it assumes that electrons are in thermal equilibrium, is 
density far from the grid (where we take = 0).] 

Ion Density rii = rioo . (1-21) 

[Applies far from grid by quasineutrality; we just assume, for the sake of this illustrative 
calculation that ion density is not perturbed by 0-perturbation.] 

Substitute: 

d 2 (j) erioc 
dx 2 eo 

This is a nasty nonlinear equation, but far from the grid \e<p/T e \ « 1 so we can use a Taylor 
expression: exp ~ 1 + ^. So 



exp (k _1 



(1.22) 



d^ en^^ eh^ 

dx 2 e Tj e Tj 1 ; 
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Solutions: <fi = (f> exp(— \x\/Xd) where 



6 ^oo 



This is called the Debye Length 

Perturbations to the charge density and potential in a plasma tend to fall off with charac- 
teristic length A/> 

In Fusion plasmas A^ is typically small, [e.g. n e = 10 20 m~ 3 T e = IkeV Ad = 2 x 10 _5 m = 
20/im] 

Usually we include as part of the definition of a plasma that A^> << the size of plasma. This 
ensures that collective effects, quasi-neutrality etc. are important. Otherwise they probably 
aren't. 



1.2.4 Plasma-Solid Boundaries (Elementary) 

When a plasma is in contact with a solid, the solid acts as a "sink" draining away the plasma. 
Recombination of electrons and ions occur at surface. Then: 

1. Plasma is normally charged positively with respect to the solid. 



Plasma 



■ 



)lid 



Sheath 

\ 



4> 



Potential 



Figure 1.7: Plasma-Solid interface: Sheath 



2. There is a relatively thin region called the "sheath", at the boundary of the plasma, 
where the main potential variation occurs. 

Reason for potential drop: 

Different velocities of electrons and ions. 

If there were no potential variation (E= 0) the electrons and ions would hit the surface at 
the random rate 

1 / 
-nv per unit area (1-25) 
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[This equation comes from elementary gas-kinetic theory. See problems if not familiar.] 



The mean speed v = 

Because of mass difference electrons move ~ faster and hence would drain out of plasma 
faster. Hence, plasma charges up enough that an electric field opposes electron escape and 
reduces total electric current to zero. 

Estimate of potential: 



Ion escape flux 
Electron escape flux 
Prime denotes values at solid surface. 
Boltzmann factor applied to electrons: 

n' e = rioo exp[e0 s /T e ] 

where <p s is solid potential relative to distant (oo) plasma. 

Since ions are being dragged out by potential assume v! i ~ 
approximately correct.] 

Hence total current density out of plasma is 



(1.26) 

1). [This is only 



1 



n'v P 



en. 



'-{vi - exp 



e<p s 



This must be zero so 



e v e 



T 1 
e In 



T e rrii 



rrii 



e 2 

[if T e = Tj.] 



For hydrogen 



e 2 

= -3.75. 

The potential of the surface relative to plasma is approximately —4 
[Note ^ is just the electron temp r in electron-volts expressed as a voltage. 



(1.27) 
(1.28) 

(1.29) 
(1.30) 
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Te 
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1.2.5 Thickness of the sheath 

Crude estimates of sheath thickness can be obtained by assuming that ion density is uniform. 
Then equation of potential is, as before, 



dx 2 



en r 



eo 



exp 




(1.31) 
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We know the rough scale-length of solutions of this equation is 



A D = (^^) 2 the Debye Length. (1.32) 

Actually our previous solution was valid only for \e<j)/T e \ << 1 which is no longer valid. 

When — e<p/T e > 1 (as will be the case in the sheath). We can practically ignore the electron 
density, in which case the solution will continue only quadratically. One might expect, 
therefore, that the sheath thickness is roughly given by an electric potential gradient 

T 1 

- -— (1.33) 
e \ D 



extending sufficient distance to reach <p s = — 4^f i.e. 

distance x ~ 4Ad 

This is correct for the typical sheath thickness but not at all rigorous. 



1.3 The 'Plasma Parameter' 

Notice that in our development of Debye shielding we used n e e as the charge density and 
supposed that it could be taken as smooth and continuous. However if the density were so 
low that there were less than approximately one electron in the Debye shielding region this 
approach would not be valid. Actually we have to address this problem in 3-d by defining 
the 'Plasma Parameter', N D , as 

N D = Number of particles in the 'Debye Sphere'. 

- //. ';tA : /, ( \ ' ' I . (1.31) 
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If N D < 1 then the individual particles cannot be treated as a smooth continuum. It will be 
seen later that this means that collisions dominate the behaviour: i.e. short range correlation 
is just as important as the long range collective effects. 

Often, therefore we add a further qualification of plasma: 

N D » 1 (Collective effects dominate over collisions) (1.35) 



1.4 Summary 

Plasma is an ionized gas in which collective effects dominate over collisions. 

[X D « size , N D » 1 .] (1.36) 
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1.5 Occurrence of Plasmas 



Gas Discharges: Fluorescent Lights, Spark gaps, arcs, welding, lighting 
Controlled Fusion 

Ionosphere: Ionized belt surrounding earth 

Interplanetary Medium: Magnetospheres of planets and starts. Solar Wind. 

Stellar Astrophysics: Stars. Pulsars. Radiation-processes. 

Ion Propulsion: Advanced space drives, etc. 

& Space Technology Interaction of Spacecraft with environment 

Gas Lasers: Plasma discharge pumped lasers: CO2, He, Ne, HCN. 

Materials Processing: Surface treatment for hardening. Crystal Growing. 

Semiconductor Processing: Ion beam doping, plasma etching & sputtering. 

Solid State Plasmas: Behavior of semiconductors. 

For a figure locating different types of plasma in the plane of density versus temperature see 
for example Goldston and Rutherford Introduction to Plasma Physics IOP Publishing, 1995, 
figure 1.3 page 9. Another is at |http: //www . plasmas . org/basics .htm 

1.6 Different Descriptions of Plasma 

1. Single Particle Approach. (Incomplete in itself). Eq. of Motion. 

2. Kinetic Theory. Boltzmann Equation. 



Odd 
dt <9x <9v 



col. 



dt 

3. Fluid Description. Moments, Velocity, Pressure, Currents, etc. 



f =%) (1-37) 



Uses of these. 

Single Particle Solutions — > Orbits 

— > Kinetic Theory Solutions — > Transport Coefs. 

— > Fluid Theory — > Macroscopic Description 
All descriptions should be consistent. Sometimes they are different ways of looking at the 
same thing. 
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1.6.1 Equations of 

V.E = 

V A E = 

F = q(E + v A B) 

1.6.2 Self Consistency 

In solving plasma problems one usually has a 'circular' system: 

The problem is solved only when we have a model in which all parts are self consistent. We 
need a 'bootstrap' procedure. 

Generally we have to do it in stages: 

• Calculate Plasma Response (to given E,B) 

• Get currents & charge densities 

• Calculate E & B for j, p. 

Then put it all together. 



Plasma Physics 

: - V.B = 

dB „ n . 1 <9E a 3 8 ) 

-at = ™ + ' ' 
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Chapter 2 

Motion of Charged Particles in Fields 



Plasmas are complicated because motions of electrons and ions are determined by the electric 
and magnetic fields but also change the fields by the currents they carry. 

For now we shall ignore the second part of the problem and assume that Fields are Prescribed. 

Even so, calculating the motion of a charged particle can be quite hard. 

Equation of motion: 

m% = q ( E + v A B ) (2.1) 

^_Jt, charge E-field velocity B-field 

Rate of change of momentum ' Lorentz Force ' 

Have to solve this differential equation, to get position r and velocity (v= r) given E(r, t), B(r, t). 
Approach: Start simple, gradually generalize. 

2.1 Uniform B field, E = 0. 

mv = qw A B (2.2) 

2.1.1 Qualitatively 

in the plane perpendicular to B: Accel, is perp to v so particle moves in a circle whose 
radius r L is such as to satisfy 

v 2 

mr L tt 2 = m— = \q\v±B (2.3) 
'7. 

Q is the angular (velocity) frequency 

1st equality shows f2 2 = v\jr\ (r L = v±/Q) 
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Figure 2.1: Circular orbit in uniform magnetic field. 



Hence second gives m^fr = \q\v±B 



m 

Particle moves in a circular orbit with 



i.e. Q = . (2.4) 



angular velocity Q = the "Cyclotron Frequency" (2-5) 

m 

v± 

and radius r { = — the "Larmor Radius. (2.6) 

2.1.2 By Vector Algebra 

• Particle Energy is constant, proof: take v. Eq. of motion then 

mv.v = — (^-mv 2 ^j = qv.(v A B) = 0. (2.7) 

• Parallel and Perpendicular motions separate. v\\ = constant because accel (oc v A B) 
is perpendicular to B. 

Perpendicular Dynamics: 

Take B in z direction and write components 

mv x = qv y B , mv y = —qv x B (2.8) 

Hence 



qB . (qB 



2 



v x = — i)y = - — v x — —Q v x (2.9) 
m \m ' 



Solution: v x = v± cos Qt (choose zero of time) 
Substitute back: 



ttl q 

Vy = -j^Vx = v±sinQt (2-10) 
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Integrate: 



x = 



v±_ . „ q v± 

s + — smSH , y = Ho ~\~ — r cos j it 



(2.11) 




q positive (ion) 
Left-handed 
Anticlockwise 

q negative (e) 
Right-handed 
Clockwise 



Figure 2.2: Gyro center (xo,yo) and orbit 

This is the equation of a circle with center r = (x , y ) and radius r L = v±/Q: Gyro Radius. 
[Angle is 6 = fit] 

Direction of rotation is as indicated opposite for opposite sign of charge: 

Ions rotate anticlockwise. Electrons clockwise about the magnetic field. 

The current carried by the plasma always is in such a direction as to reduce the magnetic 
field. 

This is the property of a magnetic material which is "Diagmagnetic" . 
When v\\ is non-zero the total motion is along a helix. 



2.2 Uniform B and non-zero E 

mv = q(E + v A B) 
Parallel motion: Before, when E = this was v\\ = const. Now it is clearly 

v... - qE W 



m 



(2.12) 



(2.13) 



Constant acceleration along the field. 



Perpendicular Motion 
Qualitatively: 

Speed of positive particle is greater at top than bottom so radius of curvature is greater. 
Result is that guiding center moves perpendicular to both E and B. It 'drifts' across the 
field. 

Algebraically: It is clear that if we can find a constant velocity that satisfies 



E + v d A B = 



(2.14) 
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Figure 2.3: E A B drift orbit 



then the sum of this drift velocity plus the velocity 

d , 



which we calculated for the E = gyration will satisfy the equation of motion. 
Take AB the above equation: 

= E A B + (v<j AB)AB = EAB + (u d .B)B - B 2 v d 



so that 

does satisfy it. 

Hence the full solution is 



E AB 

B 2 



where 



and 



V = V|| + v d + V L 

parallel cross-field drift Gyration 



1 E \\ 



rn 



(2.15) 



(2.16) 
(2.17) 



(2.18) 



(2.19) 



\d (eq |2.17 ) is the "E x B drift" of the gyrocenter. 
Comments on E x B drift: 

1. It is independent of the properties of the drifting particle (q, m, v, whatever). 

2. Hence it is in the same direction for electrons and ions. 

3. Underlying physics for this is that in the frame moving at the E x B drift E = 0. We 
have 'transformed away' the electric field. 

4. Formula given above is exact except for the fact that relativistic effects have been 
ignored. They would be important if Vd ~ c. 
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2.2.1 Drift due to Gravity or other Forces 

Suppose particle is subject to some other force, such as gravity. Write it F so that 

mv = F + gvAB = g(VvAB) (2.20) 

This is just like the Electric field case except with F/q replacing E. 
The drift is therefore 

1 F A B 

v^-p- (2.21) 

In this case, if force on electrons and ions is same, they drift in opposite directions. 

This general formula can be used to get the drift velocity in some other cases of interest (see 
later) . 

2.3 Non-Uniform B Field 

If B-lines are straight but the magnitude of B varies in space we get orbits that look quali- 
tatively similar to the E 1 B case: 



VB 




Figure 2.4: V-B drift orbit 

Curvature of orbit is greater where B is greater causing loop to be small on that side. Result 
is a drift perpendicular to both B and VB. Notice, though, that electrons and ions go in 
opposite directions (unlike E A B). 

Algebra 

We try to find a decomposition of the velocity as before into v = + v L where v rf is 
constant. 

We shall find that this can be done only approximately. Also we must have a simple expres- 
sion for B. This we get by assuming that the Larmor radius is much smaller than the scale 
length of B variation i.e., 

r L «B/\VB\ (2.22) 
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in which case we can express the field approximately as the first two terms in a Taylor 
expression: 

B ~ B + (r.V)B (2.23) 
Then substituting the decomposed velocity we get: 

= mv t = ?(vAB) = g[v L A B + v d A B + (v L + v rf ) A (r.V)B] (2.24) 
or = v d AB + v L A (r.V)B + v d A (r.V)B (2.25) 

Now we shall find that Vd/vi is also small, like r\VB\/B. Therefore the last term here is 
second order but the first two are first order. So we drop the last term. 

Now the awkward part is that Vl and r L are periodic. Substitute for r = r + r L so 

= w d A B + v L A (r L .V)B + v L A (r .V)B (2.26) 
We now average over a cyclotron period. The last term is oc e~ lQ,t so it averages to zero: 

= v rf A B + (v L A (r L .V)B) . (2.27) 

To perform the average use 



r L = (xL,y L ) = ^ (^mfit^cosfitj (2.28) 
v L = (x L ,y L ) = v± (cos fit, sin fit j (2.29) 



n 


^sinfit, 


\q\ 


v± 


^cos fit, 


\q\ 


v y y 


dy B 





So [vl A (r.V)B] s = v y y—B (2.30) 



[v L A(r.V)B] y = -v x y^B (2.31) 



(Taking VB to be in the y-direction). 
Then 



v 2 



So 



Substitute in: 



( Vy y) = - (cos fit sin fit) -± = (2.32) 

2 j 2 

(v x y) = (cos fit cos fit) ^ = 2^y^j ( 2 - 33 ) 



(v L A(r.V)B) = -j|^VS (2.34) 



= v d AB-±£vB (2.35) 
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and solve as before to get 



^VB AB „ „_> 



Vd ~ W -\q\m^^' (2 ' 36) 



or equivalently 

_ lm^ BAVB 

~~ g~2£? 52 (2 ' 37) 

This is called the 'Grad B drift'. 



2.4 Curvature Drift 

When the B-field lines are curved and the particle has a velocity v\\ along the field, another 
drift occurs. 




Center of Curvature 

Figure 2.5: Curvature and Centrifugal Force 

Take \B\ constant; radius of curvature R e . 

To 1st order the particle just spirals along the field. 

In the frame of the guiding center a force appears because the plasma is rotating about the 
center of curvature. 

This centrifugal force is F c f 

II 

F c f = 771 -zj- pointing outward (2.38) 



as a vector 



F cf = mv lj£ ( 2 - 39 ) 



[There is also a coriolis force 2m(w A v) but this averages to zero over a gyroperiod.] 
Use the previous formula for a force 

_ 1 F cf A B mvj R e A B 
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This is the "Curvature Drift". 

It is often convenient to have this expressed in terms of the field gradients. So we relate R c 
to V.B etc. as follows: 




Figure 2.6: Differential expression of curvature 



(Carets denote unit vectors) 
From the diagram 

and 

So 

But (by definition) 



db = b 2 - bi = -R c tt (2.41) 
di = oc R c (2.42) 



dX) R c R c 

~dl = ~R~ C = ~B? C 



dl 

So the curvature drift can be written 

,2 T, t-> ™„,2 



2.4.1 Vacuum Fields 



Relation between VB & R r drifts 



(2.43) 



' /b (B.V)b (2.44) 



_ m^R c B mt, B A (b,V)b 

Vd "^if A ^-^ ^ (2 ' 45) 



The curvature and V-B are related because of Maxwell's equations, their relation depends 
on the current density j. A particular case of interest is j = 0: vacuum fields. 



B w 



r{= R') 

Figure 2.7: Local polar coordinates in a vacuum field 

V A B = (static case) (2.46) 
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Consider the z-component 



1 d 

= (V A B) g = - — {rB e ) (B r = by choice). (2.47) 
r or 

° D " + * (2.48) 



dr r 

or, in other words, 



(VB) r = ~ (2.49) 



[Note also = (V A B) e = dB e /dz : (VB) Z = 0] 
and hence (V-B) pcrp = — B ^R C /R 2 C . 
Thus the grad B drift can be written: 



_ mv\ B A VB _ mVl R c A B 
VvB "^g W~~ 2q R 2 C B 2 ^ 



and the total drift across a vacuum field becomes 



1 ( 2 1 2 A R c A B 



Vr + v V b = - (m^, + -mv ± j R2B2 . (2.51) 

Notice the following: 

1. R c & Vi? drifts are in the same direction. 

2. They are in opposite directions for opposite charges. 

3. They are proportional to particle energies 

4. Curvature <-> Parallel Energy (x 2) 

V-B <-> Perpendicular Energy 

5. As a result one can very quickly calculate the average drift over a thermal distribution 
of particles because 

{\mv\) = T - (2.52) 
(-mvD = T 2 degrees of freedom (2.53) 



Therefore 



"2 



, 2TR C AB / 2TB A (b.Vjb 

(v fi + v TO ^-— U- ^ 1 (2.54) 
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2.5 Interlude: Toroidal Confinement of Single Parti- 
cles 



Since particles can move freely along a magnetic field even if not across it, we cannot ob- 
viously confine the particles in a straight magnetic field. Obvious idea: bend the field lines 
into circles so that they have no ends. 




Figure 2.8: Toroidal field geometry 



Problem 

Curvature & VB drifts 



If 2 1 2 \ RAB 

1 ( 2 1 2 \ 1 



(2.55) 
(2.56) 



Ions drift up. Electrons down. There is no confinement. When there is finite density things 




EAB 



Figure 2.9: Charge separation due to vertical drift 
are even worse because charge separation occurs — > E — > E A B — > Outward Motion. 



2.5.1 How to solve this problem? 

Consider a beam of electrons ^ fj_ = 0. Drift is 

.2 T 



Vd = 



q B T R 
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(2.57) 



What B z is required to cancel this? 

Adding B z gives a compensating vertical velocity 

B 7 



We want total 



v = uti — — for B z << B T 

Bj< 



B z mv\ q 

v z = = v\\ 1 

Bt q Bj'R 



(2.58) 



(2.59) 



So B z = —mv\\/Rq is the right amount of field. 
Note that this is such as to make 



r L {B z 



\mvu | 



R 



(2.60) 



But B z required depends on v\\ and q so we can't compensate for all particles simultaneously. 
Vertical field along cannot do it. 

2.5.2 The Solution: Rotational Transform 




Figure 2.10: Tokamak field lines with rotational transform 
Toroidal Coordinate system (r, 9, (p) (minor radius, poloidal angle, toroidal angle), see figure 



Suppose we have a poloidal field Bg 



Field Lines become helical and wind around the torus: figure 2.10 
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In the poloidal cross-section the field describes a circle as it goes round in 4>. 
Equation of motion of a particle exactly following the field is: 

Bn 



d9 B d Bq Bm 



B 



and 



r = constant. 

Now add on to this motion the cross field drift in the z direction. 




Figure 2.11: Components of velocity 



dO 

di 
dr 

dt 



Be 



v\\ + Vd cos 9 



B 

Vd sin 9 



Take ratio, to eliminate time: 



1 dr Ud sin 9 

r d9 %Lv\\ + v d cos 9 

Take Be, B,v\\,Vd to be constants, then we can integrate this orbit equation: 

[lnr] = [-ln|^i + t; d cosfl|] . 



Take r = r when cos^ = (9 = |) then 



r = r / 



i , Bvd a 
1 + cos 9 

Mil 



If J^j- << 1 this is approximately 



r = r — A cos 9 



r 



where A = J^- 
This is approximately a circular orbit shifted by a distance A: 
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Figure 2.12: Shifted, approximately circular orbit 



Substitute for v d 



B 1 mt) + \mv\) 1 

r o^-- — 2.69 

B e q v\\ B^R V ; 

1 mv\ + \mvlr p 



qB e v\\ R 

jf n A mv \\ r ° r ° 71 ^ 

UV± = ° A = WeR =rL6 R> (2 - 71) 

where r^g is the Larmor Radius in a field Bg x r/R. 

Provided A is small, particles will be confined. Obviously the important thing is the poloidal 
rotation of the field lines: Rotational Transform. 

Rotational Transform 

poloidal angle , n „. 

rotational transform = — (2.72) 

1 toroidal rotation 

(transform/27r =) t= P°^°^ a | an g| e ^ (2.73) 

toroidal angle 

(Originally, t was used to denote the transform. Since about 1990 it has been used to denote 
the transform divided by 2n which is the inverse of the safety factor.) 

'Safety Factor' 

, , 1 toroidal angle , _ , 

V a = - = — r^n — V • 2 - 74 

i poloidal angle 

Actually the value of these ratios may vary as one moves around the magnetic field. Definition 
strictly requires one should take the limit of a large no. of rotations. 

q s is a topological number: number of rotations the long way per rotation the short way. 

Cylindrical approx.: 
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In terms of safety factor the orbit shift can be written 



I A I r B ^ T 

|A| = r L6 — = r L4> —— = r L q s 



(2- 



(assuming B^ » B e ) 



2.6 The Mirror Effect of Parallel Field Gradients: E 

0, VB II B 




Figure 2.13: Basis of parallel mirror force 



In the above situation there is a net force along B. 
Force is 



< F|| > = — |gvAB| sina;= — \q\v±B sin a 
-B r 



sin a = 



B 



Calculate B r as function of B z from V.B = 0. 



Hence 



1 d d 
V.B = - — (rB r ) + —B z = 0. 
r or oz 



rB r = — I r—^-dr 



dz 



Suppose 7"l is small enough that ~ const 

[rB r ]? ~ [ 
Jo 



rL , dB z 
rdr — — 
oz 



1 2 dB z 

2 T ' L dz 



So 



r> / \ 1 dB z 



B r 



r L 1 dB z 



since = — — = +- 

B 2 2 dz 



(2- 
(2- 

(2- 
(2- 

(2- 

(2- 
(2- 
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Hence 

^ p ^_ \ n \V±r L dB g _ \mv\ dB z 

As particle enters increasing field region it experiences a net parallel retarding force. 
Define Magnetic Moment 

/I = \mv\jB . (2.85) 
Note this is consistent with loop current definition 

AT 2 l?K WlV±_ . , 

ii — AI — nr T . = 2.86 

Force is F\\ = //.VyB 

This is force on a 'magnetic dipole' of moment fj,. 

F|| = //.V||B (2.87) 

Our always points along B but in opposite direction. 

2.6.1 Force on an Elementary Magnetic Moment Circuit 

Consider a plane rectangular circuit carrying current I having elementary area alxaly = dA. 
Regard this as a vector pointing in the z direction dA. The force on this circuit in a field 
B(r) is F such that 

OB 

F x = Idy[B z (x + dx) - B z (x)\ = Idydx—^- (2.88) 

OB 

F y = -Idx[B z (y + dy)-B z (y)]=Idydx-^- (2.89) 
F z = -Idx[B y (y + dy)-B y (y)]-Idy[B x (x + dx)-B x (x)} (2.90) 
= —Idxdy 



dB x dB v 



dx dy 



Idydx—^- (2.91) 

az 



(Using V.B = 0). 

Hence, summarizing: F = IdydxV B z . Now define \i = IdA = Idydxz and take it constant. 
Then clearly the force can be written 

F = V(B.fi) [Strictly = (VB)./x] (2.92) 
H is the (vector) magnetic moment of the circuit. 

The shape of the circuit does not matter since any circuit can be considered to be composed 
of the sum of many rectangular circuits. So in general 

/i = IdA (2.93) 
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and force is 

F = V(B./i) (fi constant), (2.94) 

We shall show in a moment that is constant for a circulating particle, regard as an 
elementary circuit. Also, fi for a particle always points in the -B direction. [Note that this 
means that the effect of particles on the field is to decrease it.] Hence the force may be 
written 

F = -fi\7B (2.95) 



This gives us both: 

• Magnetic Mirror Force: 

and 

• Grad B Drift: 



F\\ = -nV\\B (2.96) 



1 FAB /iBAVfl , on _, 



2.6.2 fi is a constant of the motion 

'Adiabatic Invariant' 



Proof from Fu 

Parallel equation of motion 



dv\\ dB 



So 



or 



dvu d ,1 2 \ dB dB 

mv\\—r- = -r(-mvu) = —uv z —— = —u—r- (2.99) 
11 dt dt K 2 11 ; p dz p dt K ' 

d . 1 on dB 



Conservation of Total KE 



Combine 



Y{\ mv l + \ mv l) = (2.101) 

(XL Zi Zi 

^-mv\ + siB) = a (2.102) 



= ^ = As required (2.104) 
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Angular Momentum 



of particle about the guiding center is 

mv± 2m \mv\ 

r L mv± = T - rF rmv± = — r „ (2.105) 
\q\B \q\ B 

= ^/i . (2.106) 

m 

Conservation of magnetic moment is basically conservation of angular momentum about the 
guiding center. 

Proof direct from Angular Momentum 

Consider angular momentum about G.C. Because 6 is ignorable (locally) Canonical angular 
momentum is conserved. 

p = [r A (mv + qA)] z conserved. (2.107) 

Here A is the vector potential such that B = V A A 
the definition of the vector potential means that 

ft = 1 (2.08) 
r or 



T 2 

r L A e (r L ) = [ L r.B z dr= r ^B z = ^ (2.109) 
Jo 2 \q\ 



Hence 



— q mfi 

p = -j— 7-rLV±m + q— — (2.110) 



in \q 
\q 



-T^m/x. (2.111) 



So p = const = constant. 

Conservation of /x is basically conservation of angular momentum of particle about G.C. 
2.6.3 Mirror Trapping 

F|| may be enough to reflect particles back. But may not! 
Let's calculate whether it will: 
Suppose reflection occurs. 
At reflection point v\\ r = 0. 
Energy conservation 

2™( u lo + 4>) = 2 mu lr (2-112) 
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Reflection 




B-lines 



Stronger 



pi conservation 



Hence 



Figure 2.14: Magnetic Mirror 



k mv ±o 



\mv\ r 



Bo 



B r 



2 , 2 
V ±0 + V \\0 = 


B r 2 
Bo 


Bo 




B r 





2.6.4 Pitch Angle 9 



tan 9 

Bo 

B, r 



v± 

V\\ 



v lo + v \\o 



sin 2 6*o 



So, given a pitch angle 9q, reflection takes place where B /B r = sin 2 9q. 
If 6q is too small no reflection can occur. 
Critical angle 9 C is obviously 

9 C = sin' 1 (B /B 1 )^ 

Loss Cone is all 9 < 9 C . 

Importance of Mirror Ratio: R m = Bi/B . 



2.6.5 Other Features of Mirror Motions 

Flux enclosed by gyro orbit is constant. 



2 2 

■Km v , 



$ = uriB = 

q 2 B 2 
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B 



Figure 2.15: Critical angle 9 C divides velocity space into a loss-cone and a region of mirror- 
trapping 



9 

Note that if B changes 'suddenly' ji might not be conserved 



MS 



Flux Tube 




(2.120) 



2-nm \mv\ 
2nm 

— —[/, = constant. (2.121) 



Figure 2.16: Flux tube described by orbit 

Basic requirement 

r L « B/\VB\ (2.122) 

Slow variation of B (relative to r L ). 

2.7 Time Varying B Field (E inductive) 

Particle can gain energy from the inductive E field 

<9B 

V A E = -— (2.123) 
or j E.dl = -^B.ds = -^ (2.124) 

Hence work done on particle in 1 revolution is 

/( • d<& 
\q\E.d£ = +\q\ J B.ds = +\<l\^ = \q\Bmj, (2.125) 
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Figure 2.17: Particle orbits round B so as to perform a line integral of the Electric field 
(d£ and vj_g are in opposition directions). 

= 1^-1 = ^^ (2-126) 



Hence 



but also 



Hence 



d (\ 9 \ M /l ,\ cfc 

sM-krHa™ 4 )""* (2 - 128) 



|(^) = |(MB). (2.29) 
| = 0. (2.130) 



Notice that since $ = ^pr/i, this is just another way of saying that the flux through the gyro 
orbit is conserved. 

Notice also energy increase. Method of 'heating'. Adiabatic Compression. 



2.8 Time Varying E-field (E, B uniform) 

Recall the E A B drift: 



EAB 

when E varies so does veab- Thus the guiding centre experiences an acceleration 



VEAB = (2.131) 



d /EAB\ 



In the frame of the guiding centre which is accelerating, a force is felt. 

d /EAB 

dt I - B 



F a = —m— ( — — — J (Pushed back into seat! — ve.) (2.133) 
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This force produces another drift 



1 F a A B m d /EAB\ 

= - — ^f(E.B)B- J B 2 E) 

qB dt^ J > 

m ■ 



qB 2 



This is called the 'polarization drift 7 . 



(2.134) 

(2.135) 
(2.136) 



Vfl = VE/\B + Vp = 



E AB 



m 



B 2 qB 2 
E AB 1 



B 2 



QB 



+ — E 
E 



(2.137) 
(2.138) 




Initiar^osition at rest (when E = 0). Final average (gyro) locus. 



Figure 2.18: Suddenly turning on an electric field causes a shift of the gyrocenter in the 
direction of force. This is the polarization drift. 

Start-up effect: When we 'switch on' an electric field the average position (gyro center) of 
an initially stationary particle shifts over by ~ \ the orbit size. The polarization drift is this 
polarization effect on the medium. 

Total shift due to v p is 



Ar / Vpdt 



m 



qB 2 



E i dt 



m 



qB 2 



[AE^ 



(2.139) 



2.8.1 Direct Derivation of ^ effect: 'Polarization Drift 



dt 



Consider an oscillatory field E = Ee twt 

dv 



m 



dt 



(± r B) 

q (E + v A B) 
q (Ee" 



+ v A 



Try for a solution in the form 



v = v D e-^' + v L 



(2.140) 
(2.141) 

(2.142) 
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where, as usual, v L satisfies mw L = q\ L A B 
Then 



-iuit 



(1) m(-icuv D = q (E + v D A B) xt 
Solve for v D : Take AB this equation: 

(2) - miuj (v D A B) = q (E A B + (B 2 .v|£>) B - B 2 v D 
add miuj x (f) to q x (2) to eliminate v D A B. 

m 2 uj\ D + q 2 (E A B - £ 2 vd) = miuoqE 

E AB 



or : 



i.e. 



Vfl 



f - 



2 2 

q 2 B 2 

n 2 



miu 
iuiq 



E + 



fi 2 
E AB 



(2.143) 
(2.144) 

(2.145) 
(2.146) 

(2.147) 



Since — iu 



dt 



this is the same formula as we had before: the sum of polarization and 



E A B drifts except for the [1 — uj 2 VL 2 } term. 
This term comes from the change in with time (accel). 
Thus our earlier expression was only approximate. A good approx if uj « Q. 



2.9 Non Uniform E (Finite Larmor Radius) 



m 



dt 



q (E(r) + v A B) 



Seek the usual soltuion v = v D + v 9 . 
Then average out over a gyro orbit 

dv 



Hence drift is obviously 



(m^) = = (q (E(r) + v A B)) 
= g[(E(r))+v D AB] 

(E(r))AB 
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So we just need to find the average E field experienced. 
Expand E as a Taylor series about the G.C. 

/ x 2 d 2 v 2 d 2 \ 
E(r) = E + (r.V) E + ( — + | — J E + cross terms + 



(2.148) 



(2.149) 
(2.150) 

(2.151) 



(2.152) 
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(E.g. cross terms are xy^ y 'E). 
Average over a gyro orbit: r = r L (cos9, sin6*, 0). 
Average of cross terms = 0. 
Then 



(E(r))=E + ((r L >.V)E + ^V 2 E. ( 



linear term (n) = 0. So 



2 

(E(r)) ~ E + T -±V 2 E ( 



Hence E A B with 1st finite-Larmor-radius correction is 

/ r| „\EAB 
^AB=(^l + fV 2 J^^. ( 

[Note: Grad B drift is a finite Larmor effect already.] 
Second and Third Adiabatic Invariants 

There are additional approximately conserved quantities like \i in some geometries. 



2.10 Summary of Drifts 



E A B 

v £ — — Electric Field 
B 2 

1FAB n IP 

v F = — — General Force 

q B 2 

(j.2 \ E A B 
1 + — V 2 Nonuniform E 
4 ) B 2 

mv\ B A VB _ _ 

vvb = -z ™ — GradB 

2q B 6 

mvf R c A B 
vr = „ 9 „ 9 Curvature 

1 / 2 1 2 \R C AB „ , , 

Vfl + vvB = - rra;,, + -mu 1 n9 Vacuum Fields, 

g V 11 2 / /i 2 !* 2 

o E_\_ 

v p = — Polarization 



\q\ \n\B 



Mirror Motion 



Force is F = —jiVB. 



\i = is constant ( 

2B 
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Chapter 3 



Collisions in Plasmas 



3.1 Binary collisions between charged particles 

Reduced-mass for binary collisions: 

Two particles interacting with each other have forces 

F12 force on 1 from 2. 

F 2 i force on 2 from 1. 

By Newton's 3rd law, F 12 = — F 21 . 

Equations of motion: 



Combine to get 
which may be written 



mifi = F12 ; m 2 r 2 = F 21 (3.1) 



fx - r 2 = F 12 (— + — ) (3.2) 
\m 1 m 2 



mim 2 d 2 

(ri - r 2 ) = F12 (3.3) 



mi + m 2 dt 2 

If F 12 depends only on the difference vector ri — r 2 , then this equation is identical to the 
equation of a particle of "Reduced Mass" m r = moving at position r = r x — r 2 with 

respect to a fixed center of force: 

m r r = Fi 2 (r) . (3.4) 

This is the equation we analyse, but actually particle 2 does move. And we need to recognize 
that when interpreting mathematics. 

If F21 and ri — r 2 are always parallel, then a general form of the trajectory can be written 
as an integral. To save time we specialize immediately to the Coulomb force 

F U = P 4 (3.5) 

Solution of this standard (Newton's) problem: 
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1: projectile 



b Impact Parameter' 

2: target 

Figure 3.1: Geometry of the collision orbit 



Angular momentum is conserved: 

m r r 2 9 = const. = m r bv\ 

Substitute u=\ then Q= b -^ = u 2 bv x 
Also 

d 1 



clockwise from symmetry) 



1 du 



dt u u 2 d9 

d 2 u ■ 



-bvi 



du 

de 



-bv% 



-{bvxf u 2 



d6 2 



<fu 

de 2 



Then radial acceleration is 



f — r6 2 = — (bv i) 2 u 2 



'dtu 



+ u 



'12 



i.e. 



d 2 u 



+ u 



dO 2 ' " 4vre m r {bvif 
This orbit equation has the elementary solution 



_ 1 n a ™ 2 
u = - = C cos — 



1 



47re m r (bvi)' 



(3.6) 



(3.7) 
(3.8) 

(3.9) 
(3.10) 

(3.11) 



The sin# term is absent by symmetry. The other constant of integration, C, must be deter- 
mined by initial condition. At initial (far distant) angle, 9i, u\ — — = 0. So 



= Ccos9 1 



Q1Q2 



47re m r {bvi) A 



There: 



Hence 



ri = —v\ = —ovi — U = +bviC sm&i 
dd 



tan 6*i 



-1/Cb 



sin 6\ 

1 Aneo mr(bvi) 1 



b 



(3.12) 

(3.13) 
(3.14) 
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where 

»90 = -; 2 • ( 3 - 15 ) 

47re m r vf 

Notice that tan0 1 = — 1 when 6 = 6 90 . This is when 6\ = —45° and x — 90°. So particle 
emerges at 90° to initial direction when 

b = b 90 "90° impact parameter" (3.16) 

Finally: 

1 1 / h 2 \ 5 

C =-- b COSeCdl = --b{ 1 + f) (3 ' 17) 

3.1.1 Frames of Reference 

Key quantity we want is the scattering angle but we need to be careful about reference 
frames. 

Most "natural" frame of ref is "Center-of-Mass" frame, in which C of M is stationary. C of 
M has position: 

R= m iri + m 2 r 2 
mi + m 2 



and velocity (in lab frame) 



Now 



v= m lVl+ m 2 v 2 
mi + m 2 



n = R+ m2 r (3.20) 
mi + m 2 

r 2 = R r (3.21) 

mi + m 2 

So motion of either particle in C of M frame is a factor times difference vector, r. 
Velocity in lab frame is obtained by adding V to the C of M velocity, e.g. m ^^ 2 + V. 
Angles of position vectors and velocity differences are same in all frames. 
Angles (i.e. directions) of velocities are not same. 

3.1.2 Scattering Angle 

In C of M frame is just the final angle of r. 

-26 l + X = n (3.22) 

(9i is negative) 

X = n + 29 1 ; 6 1 = . (3.23) 
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Figure 3.2: Relation between 9\ and \. 



tan 0! = tan = ^-\ = - cot | (3.24) 

So 

cot I = A (3.25) 

tan | = ^ (3.26) 

But scattering angle (defined as exit velocity angle relative to initial velocity) in lab frame 
is different. 

Final velocity in CM frame 

v cm = wicm (cosXc,sinxc) = — ^ — v i ( cos Xc sin Xc) (3.27) 

mi + m 2 

[ Xc = X an d v i is initial relative velocity]. Final velocity in Lab frame 

< = v'cm + V = (V + cos Xc , sin Xe ) (3.28) 
V mi + m 2 mi+m 2 J 

So angle is given by 

v + ^rk cos Xc v mi + m 2 

cot XL = 1+ ? = cosecx c + cot x c 3.29 

, 1 sin y t f i m 2 

For the specific case when m 2 is initially a stationary target in lab frame, then 

V = mit>1 — and hence (3.30) 
mi + m 2 

cotxi = — -cosecXc + cot Xc (3.31) 
m 2 

This is exact. 

Small angle approximation (cot% — > ^,cosecx — > ^ gives 

1 mi 1 1 m 2 

— = —- + - ^ Xl = 3.32 

Xl m 2 Xc Xc m 1 + m 2 

So small angles are proportional, with ratio set by the mass-ratio of particles. 
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Figure 3.3: Collisions viewed in Center of Mass and Laboratory frame. 



3.2 Differential Cross-Section for Scattering by Angle 

Rutherford Cross-Section 

By definition the cross-section, a, for any specified collision process when a particle is passing 
through a density n 2 of targets is such that the number of such collisions per unit path length 
is n 2 a. 

Sometimes a continuum of types of collision is considered, e.g. we consider collisions at 
different angles (x) to be distinct. In that case we usually discuss differential cross-sections 
(e.g ^) defined such that number of collisions in an (angle) element d\ per unit path length 
is n 2 ^dx. [Note that ^ is just notation for a number. Some authors just write but I 

find that less clear.] 

Normally, for scattering-angle discrimination we discuss the differential cross-section per unit 
solid angle: 

do 

This is related to scattering angle integrated over all azimuthal directions of scattering by: 
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Figure 3.4: Scattering angle and solid angle relationship. 



dQ s = 27T sin x<ix (3.34) 

So that since 

da ^ da , , 

dQ s = —dx (3.35) 



dCl s dx 
we have 

da 1 da , 

f Q Oft \ 

dfl s 27rsinx dx 

Now, since x is a function (only) of the impact parameter, b, we just have to determine the 
number of collisions per unit length at impact parameter b. 




Figure 3.5: Annular volume corresponding to db. 



Think of the projectile as dragging along an annulus of radius b and thickness db for an 
elementary distance along its path, d£. It thereby drags through a volume: 

d£27cbdb . (3.37) 

Therefore in this distance it has encountered a total number of targets 

d£2nbdb . n 2 (3.38) 
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at impact parameter b(db). By definition this is equal to d£^dbn 2 . Hence the differential 
cross-section for scattering (encounter) at impact parameter b is 

2tt6 . 



da 
db 



Again by definition, since x is a function of b 



da 
dx 



dx 



da 
db' 



db 



da 
dx 



da 

~db 



db 
dx 



[db/dx is negative but differential cross-sections are positive.] 
Substitute and we get 



da 
dQc 



1 da 
2n sin x db 



db 
dx 



b 

sin 



db 
dx 



[This is a general result for classical collisions.] 
For Coulomb collisions, in C of M frame, 



db 
dx 



= b 



90 



cot 

d 

dx 



b 

bgo 



, X 
cot — = 

2 



090 2 X 

— cosec — 
2 2 



Hence 



da 
dtt s 



30 cot f b 90 



2 X 

— cosec — 
2 2 



sin \ 

2 
^90 

2 2 sin | cos | sin 2 | 



bin cos |/ sin | 



bio 



4 sin 4 f 



This is the Rutherford Cross-Section. 



da 
dQ* 



b 2 



4 sin 4 f 



(3.39) 



(3.40) 



(3.41) 



(3.42) 
(3.43) 

(3.44) 
(3.45) 
(3.46) 



(3.47) 



for scattering by Coulomb forces through an angle x measured in C of M frame. 
Notice that Jf — > oo as x 0. 

This is because of the long-range nature of the Coulomb force. Distant collisions tend to 
dominate, (x — > <^ b — > oo). 
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3.3 Relaxation Processes 



There are 2 (main) different types of collisional relaxation process we need to discuss for a 
test particle moving through a background of scatterers: 

1. Energy Loss (or equilibrium) 

2. Momentum Loss (or angular scattering) 

The distinction may be illustrated by a large angle (90°) scatter from a heavy (stationary) 
target. 

If the target is fixed, no energy is transferred to it. So the energy loss is zero (or small if 
scatterer is just 'heavy'). However, the momentum in the x direction is completely 'lost' in 
this 90° scatter. 

This shows that the timescales for Energy loss and momentum loss may be very different. 
3.3.1 Energy Loss 

For an initially stationary target, the final velocity in lab frame of the projectile is 

, / mvi m 2 v 1 m 2 v 1 . \ 

Vt = 1 cos Xc , sm Xc (3.48) 

\m\ + m 2 m 1 + m 2 m 1 + m 2 J 

So the final kinetic energy is 

1 /2 1 2 \ ( mi \ 2 2m\m 2 



K' = miv'l = - mi v 2 J +- cosxc (3.49) 

2 2 {\m 1 + m 2 J ( mi + m 2 ) 



+ 7 (cos 2 x c + sin 2 Xc ) \ (3.50) 

(mi + m 2 J v 7 J 

= \™<ivl\l+ 2m i m 2 (cosxc-m (3.51) 
2 { (mi + m 2 ) J 



1 m .,2 Jl , J^o„ ;J Xc 

-mi^i < 1 + ■ ^ 

* { (mi + m 2 



2 I ^/;ti/;t 2 . 2 Ac I / Q ko x 

mi^i < 1 H 2 2 sin — > (3.52) 



Hence the kinetic energy lost is AK = K — K' 

1 4mim 2 . 2 Xc f0 „v 

= ^if 2 t ; T2 sm — (3.53) 

* (mi + m 2 J z 

1 2 4mim 2 1 r . ,Xcb 

= -miV 1 2 2 [using cot — = - — J (3.54) 

2 (mi + m 2 ) ( A-) _i_ i 2 o 90 

(exact). For small angles x ^ 1 i- e - b/b$ ^> 1 this energy lost in a single collision is 
approximately 



1 2 \ 4mim 2 /&9o 



2 



2 ^ ; I 1 \ 2 7 (3.55) 
2 / (mi + m 2 ) \ o / 
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If what we are asking is: how fast does the projectile lose energy? Then we need add up the 
effects of all collisions in an elemental length dl at all relevant impact parameters. 

The contribution from impact parameter range db at b will equal the number of targets 
encountered times AK: 

jdo ujh 1 2 4m i™2 (b 9 oY ,„ - R v 

^dzznbab -m\v\ ^ ~r (3.56) 



encounters 



2 (mi + m 2 ) \ b 



Loss per encounter (AK) 
This must be integrated over all b to get total energy loss. 

Am l m 2 f fb 90 x2 



dK = n 2 d£K '—^ / Mr 2nhdh ( 3 - 57 ) 

(mi + m 2 ) J \ b J 

so 

^ = ^ % ~~~~~T2 Stt^ [In (3.58) 
cm- (mi + m 2 ) 

We see there is a problem both limits of the integral (6 — > 0, & — > oo) diverge logarithmically. 
That is because the formulas we are integrating are approximate. 

1. We are using small-angle approx for AK. 

2. We are assuming the Coulomb force applies but this is a plasma so there is screening. 

3.3.2 Cut-offs Estimates 

1. Small- angle approx breaks down around b = 690 . Just truncate the integral there; 
ignore contributions from b < 6 90 . 

2. Debye Shielding says really the potential varies as 



exp 3— 1 

oc ^1 instead of oc - (3.59) 

r r 

so approximate this by cutting off integral at 6 = Ad equivalent to 

&min = &90- &max = ^D- (3.60) 

f = ^^PV^blAI (3.61) 
M (mi + m 2 ) 



A = -— = -— ^ ) / — (3.62) 
90 \ne z J / \47re m r t>i/ 
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So Coulomb Logarithm is 'In A' 

A 



2 



QlQ2 



690 \ne 2 J I \Ait€om r v l/ 

Because these cut-offs are in In term result is not sensitive to their exact values. 
One commonly uses Collision Frequency. Energy Loss Collision Frequency is 



Substitute for 6 90 and m T (in 690) 



1 dK m\m2 n ,0 . 1 . 1 

n 2 v 1 ■ -j 8tt6 90 In |A| 

(mi + m 2 ) 



v K 



n 2 v 1 



17111712 



m 1 + m 2 y 

8tt 



8tt 



9i 92 



4 7 re JM ^' 2 



' mi+r«2 



In A 



n2 



qIq 2 



(47re ) 2 mim 2 vl 

Collision time Tk = 1/vk 
Effective (Energy Loss) Cross-section 



In A 



1 dK 

K dl 



a K n 2 



o~k = VK/n 2 Vi = 



8tt 



(47re ) rniin2vf 



In A 



(3.63) 



(3.64) 



(3.65) 
(3.66) 



(3.67) 



3.3.3 Momentum Loss 

Loss of x-momentum in 1 collision is 



= m 1 (f 1 
= m\V\ jl - 

= Px 



J Lxi 



m 2 



m x + m 2 



+ 

mi + rn2 mi + iri2 

(1 - COSXc) 



COSXc 



)} 



Px 



m 2 



xi 



nil + 2 



Px 



m 2 



2h 2 



nil + itt-2 b 2 



(3.68) 
(3.69) 

(3.70) 
(3.71) 



(small angle approx). Hence rate of momentum loss can be obtained using an integral 
identical to the energy loss but with the above parameters: 



dp 
dl 



n 2 p 
n 2 p 



m2 



b mm ml+m 2 b 2 

™2 , ,2 



2b 2 

Zt)9 ° 2«bdb 



nil + m 2 



4n bg In A 



Note for the future reference: 

dp 



dp 2 niim 2 

"tt = v ~n> = n i v 

at di nil + rn,2 



Anbl In A. 



(3.72) 
(3.73) 

(3.74) 



49 



Therefore Momentum Loss. 
Collision Frequency 

1 dp m 2 



v " = Vl ~~o = n2Vl ~ — V~ZT ^^90 my ^ (3.75) 

1 2 



p d£ mi + m 2 



m 2 

n 2 Vi 47r 

mi + m 2 



glg2 
u mi+m2 1 -I 



In A (3.76) 



n q i q % 4:71 ( mi + m ^ l n A (377) 
" 2 (4vre ) 2 m 2 m\v\ 



Collision Time r p = l/v p 
Cross-Section (effective) a = u p /n 2 vi 
Notice ratio 



Energy Loss is K 2 /m 1 4- m 2 2mi 



Momentum loss u p m 1 m 2 j m 2 m 2 mi + m 2 

This is 



(3.78) 



~ 2 if mi » m 2 (3.79) 
= 1 if nii = m 2 (3.80) 
< 1 if mi« m 2 . (3.81) 

Third case, e.g. electrons — > shows that mostly the ang/e of velocity scatters. Therefore 
Momentum 'Scattering' time is often called '90° scattering' time to 'diffuse' through 90° in 
angle. 

3.3.4 'Random Walk' in angle 

When nil « m 2 energy loss << momentum loss. Hence \v' L ~ Vi. All that matters is the 
scattering angle: xl-Xc- 2b 90 /b. 

Mean angle of deviation in length L is zero because all directions are equally likely. 
But: 

Mean square angle is 



Aa 2 = n 2 L / x Znbdb (3.82) 
= Ln 2 8tt b 2 90 In A (3.83) 



Spread is 'all round' when Aa 2 ~ I. This is roughly when a particle has scattered 90° on 
average. It requires 

Ln 2 Snb 2 90 In A = 1 . (3.84) 



50 



So can think of a kind of 'cross-section' for £ <7 90 ' 90° scattering as such that 

n 2 L'o"9o' = 1 when Ln 2 Sn 6g In A = 1 (3.85) 
i.e. 'a 90 ' = 87r6g lnA (= 2a p ) (3.86) 
This is 8 In A larger than cross-section for 90° scattering in single collision. 

Be Careful! 'a^ is not a usual type of cross-section because the whole process is really- 
diffusive in angle. 

Actually all collision processes due to coulomb force are best treated (in a Mathematical 
way) as a diffusion in velocity space 

— > Fokker-Planck equation. 
3.3.5 Summary of different types of collision 

The Energy Loss collision frequency is to do with slowing down to rest and exchanging 
energy. It is required for calculating 

Equilibration Times (of Temperatures) 
Energy Transfer between species. 
The Momentum Loss frequency is to do with loss of directed velocity. It is required for 
calculating 

Mobility: Conductivity/Resistivity 

Viscosity 

Particle Diffusion 

Energy (Thermal) Diffusion 
Usually we distinguish between electrons and ions because of their very different mass: 

Energy Loss [Stationary Targets] Momentum Loss 



K 



K 



K 



Vi, 



n P 



57T 

\2 rr - l 2, il 3 



Vpi = n, 



(47re ) m 2 e v e 
Z 2 e 4 8tt 

(47re ) 2 m e miV : 
Z 2 e 4 8tt 



In A 



In A p u P 



K 



K 



(47re ) 2 m 2 vf 
Z 2 e A 



In A 



m e + m e 
2m e 

m e + mi 

2m e 
mi + TOj 



2m P 



1 



7T 



n P 



(47re ) mim e vf 



In A v v iP = K Vi, 



x 



2rrii 
m e + rrii 1 
2mij 2 



(3.87) 



Sometimes one distinguishes between 'transverse diffusion' of velocity and 'momentum loss'. 
The ratio of these two is 



p 2 AL. 





dxl 1 


1 dp 


pAL 


dL 1 


p dL 




( mi 


Xc) 2 




\mi+m2 




TO2 


xl - 



m\+rri2 2 



2m 2 
m 1 + m 2 ' 



(3.88) 
(3.89) 



51 



So 



Hence 



'0-90' 2m 2 



'o-p' m 1 + m 2 



= 1 like particles (3.90) 

~ 2 mi« m 2 (3.91) 

~ — m 2 « mi . (3.92) 

777,1 



= V= \e(='^'!!) (3.93) 

V* = 2 p v ei = V-2 2 (=^ £ e)(=V a ') (3.94) 

77/g 

± i/jj = p h>u = K va (= uaW) (Like Ions) (3.95) 

_L 2m e m e k K /o nr\ 

Vie = P V ie = V ie = I/jj = I/jj (3.96) 

777j 777j 

[But note: ions are slowed down by electrons long before being angle scattered.] 



3.4 Thermal Distribution Collisions 

So far we have calculated collision frequencies with stationary targets and single-velocity 
projectiles but generally we shall care about thermal (Maxwellian) distributions (or nearly 
thermal) of both species. This is harder to calculate and we shall resort to some heuristic 
calculations. 



3.4.1 



Very rare for thermal ion velocity to be ~ electron. So ignore ion motion. 
Average over electron distribution. 

Momentum loss to ions from (assumed) drifting Maxwellian electron distribution: 

\2" 



/e(v) = ne v2S;) 2 exp 



m(y- v d ) 
2T 



(3.97) 



Each electron in this distribution is losing momentum to the ions at a rate given by the 
collision frequency 

q 2 e qj 47T (m e + m, 



Hi 



(47re r 



rriimlv 3 



In A 



so total rate of loss of momentum is given by (per unit volume) 

dp 
~~ dt 

To evaluate this integral approximately we adopt the following simplifications. 



J /e(v) v p {v) m e v d\ 



(3.98) 



(3.99) 
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1. Ignore variations of In A with v and just replace a typical thermal value in A = 

Ad/&90(^i)- 

2. Suppose that drift velocity v d is small relative to the typical thermal velocity, written 
v e = \Jt e /m e and express f e in terms of u = ^ to first order in u rf e 



fe 



n,. 



(27T) 2 V 



cxp 



- 1 r ^ 



1 + u.u d ) exp 



1 -|- U x ll(ij f 



taking x-axis along and denoting by f Q the unshifted Maxwellian. 
Then momentum loss rate per unit volume 



d Px f , ,3 

= / f e v v m e v x d v 



dt 



/V 
(l + u x u d )f -| v x d 3 v 

/U 2 
fo d 3 v 



To evaluate this integral, use the spherical symmetry of f Q to see that: 



6 JO V 

/ /oS 
JO 



3 
2tt 



A^v 2 dv 
Jlvdv 



n. 



3 (27r)5^^o 
2vr n e _ 
3 (2tt)2 



cxp 



— V 

2rf 



dv 2 



2 



Thus the Maxwell-averaged momentum-loss frequency is 

1 dp 2 
p dt 



V e i = 



3(2tt 

j 



3(2tt 

(where p = m e Vdn e is the momentum per unit volume attributable to drift). 



r n i 

3(2tt)2 (47re 



mim\v\ 



In A, 



2 / ze 2 \ 2 4tt 



3(2tt) 



n 7 . 



47re n 



, — t lnA e 



(3.100) 
(3.101) 



(3.102) 



(3.103) 



(3.104) 



(3.105) 



(3.106) 



(substituting for thermal electron velocity, v e , and dropping ^ order term), where Ze = ^. 
This is the standard form of electron collision frequency. 
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3.4.2 i -> e 



Ion momentum loss to electrons can be treated by a simple Galilean transformation of the 
e —>■ i case because it is still the electron thermal motions that matter. 

Ions — ► Electrons i ons < — Electrons 



. • . . = • . • * 




u» ^ = v d = -v< 



Figure 3.6: Ion-electron collisions are equivalent to electron-ion collisions in a moving refer- 
ence frame. 



Rate of momentum transfer, is same in both cases: 

-pv (3.107) 



■' df> 

dp 
dt 

Hence p e v ei = piV ie or 



_ Pe_ n e m e _ 

V ie = —Vei = Vei (3.108) 

Pi riirrii 

(since drift velocities are the same). 

Ion momentum loss to electrons is much lower collision frequency than e — ^ i because ions 
possess so much more momentum for the same velocity. 



3.4.3 i — > i 

Ion-ion collisions can be treated somewhat like e — > i collisions except that we have to 
account for moving targets i.e. their thermal motion. 

Consider two different ion species moving relative to each other with drift velocity v^', the 
targets' thermal motion affects the momentum transfer cross-section. 

Using our previous expression for momentum transfer, we can write the average rate of 
transfer per unit volume as: [see 3.74 "note for future reference"] 

-lt= f f y r **** & 90 In A hhd 2 Vl d\ 2 (3.109) 

at J J mi + m 2 



where v r is the relative velocity (vi — v 2 ) and 690 is expressed 

Q1Q2 1 



^90 



47re m r v 2 r 
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(3.110) 



and m r is the reduced mass mi J 7 " 2 . 

' mi+m2 



Since everything in the integral apart from j\f 2 depends only on the relative velocity, we 
proceed by transforming the velocity coordinates from Vi, v 2 to being expressed in terms of 
relative (v r ) and average (V say) 



vi - v 2 



V 



TO1V1 + m 2 v 2 



mi + m 2 

Take fi and / 2 to be shifted Maxwellians in the overall C of M frame: 

\2l 



(3.1H) 



( m 3 \ 

\2tiTJ 



fj = n j 

where miV dl + m 2 v d2 = 0. Then 



cxp 



2T 



hh 



nin 2 



( mi \2 / m 2 \ 
\2ttt) \2ttt) 



cxp 



(J = l,2) 



m\v\ m 2 v\ 



(3.112) 



2T 



2T 



(3.113) 



to first order in v d . Convert CM coordinates and find (after algebra) 
hh = HlU2 (2^) 2 (2SO 2 6XP 



~~2T 2T~ 



x ^ 1 + — v d .v, 



(3.114) 



where M = m x + m 2 . Note also that (it can be shown) d 3 v 1 d 3 v 2 = d 3 v r d 3 V. Hence 



dp 

~dt 



II 



\ r m r v r A-K bl Q mAnin 2 



M \ 2 ( m 



(m r \ 



2itTJ \2nTj 



exp 



MV 2 
2T 



cxp 



m r v 2 r 
2T 



){l + yV,v r }rfVV (3.115) 



and since nothing except the exponential depends on V, that integral can be done: 
v r m r v r 4ir In A n\n 2 



dp 

~dt 



\2irTj 



cxp 



-m r v^ 



2vr 



m r 



l+'-^v d .v r \d 3 v r (3.116) 



This integral is of just the same type as for e — i collisions, i.e. 



dp 
dt 



v d v rt m r An b go (v rt ) lnA t n x n 2 \ — f {v r )d 6 v r 

U3 



v d v rt m r A n bl (v rt ) In A t n x n 2 



3(2tt) 



(3.117) 
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where v rt = y— , bc, (v r t) is the ninety degree impact parameter evaluated at velocity v tr , 
and f is the normalized Maxwellian. 

dp 2 (qiq 2 \ 2 4tt 

- — = r ——5- In A t nin 2 m r v d (3.118) 

3(27r)2 V47re / mf,v?. t 

This is the general result for momentum exchange rate between two Maxwellians drifting at 
small relative velocity 

To get a collision frequency is a matter of deciding which species is stationary and so what the 
momentum density of the moving species is. Suppose we regard 2 as targets then momentum 
density is n\m\Vd so 



1 dp 

V\2 = 



mJm 2 ^-^ . (3.H9) 

2 \47reo/ m r Vrt rri\ 



n x m x v A dt 3(2vr)2 \47re / m r vf t 

This expression works immediately for electron-ion collisions substituting m r ~ m e , recov- 
ering previous. 

2 J I 

For equal-mass ions m r = — = \vrii and v rt = J — = J—. 
Substituting, we get 

Va = A m (|^Y 4^ In A (3.120) 
37T2 \47re / m ?j"2 

that is, 4^ times the e — i expression but with ion parameters substituted. [Note, however, 
that we have considered the ion species to be different.] 



3.4.4 e -> e 

Electron-electron collisions are covered by the same formalism, so 



1 / e 2 \ 2 4tt 

V ee = -^nA- -^-r InA . (3.121) 

3tt2 V 47re o/ miTe 2 

However, the physical case under discussion is not so obvious; since electrons are indistigu- 
ishable how do we define two different "drifting maxwellian" electron populations? A more 
specific discussion would be needed to make this rigorous. 

Generally v ee ~ is ei /\/2 : electron-electron collision frequency ~ electron- ion (for momentum 
loss) . 

3.4.5 Summary of Thermal Collision Frequencies 

For momentum loss: 

"a = ^ rii (^-) — ^ lnA e . (3.122) 
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v ee ~ — ^= z/ ei . (electron parameters) (3.123) 

F ie = ^v ei . (3.124) 

Fjj, = y~ rtn ( \ — i — g- ( — ^ InAj (3.125) 

3^71" V47re / m ?y.2 Vmi+m,,/ 



Energy loss K v related to the above ( p u) by 



K v = 2m *v . (3.126) 
mi + m 2 



Transverse 'diffusion' of momentum - L i/, related to the above by: 



v = 2m2 *v . (3.127) 
mi+m 2 



3.5 Applications of Collision Analysis 

3.5.1 Energetic ('Runaway') Electrons 

Consider an energetic (|m e t>^ >> T) electron travelling through a plasma. It is slowed down 
(loses momentum) by collisions with electrons and ions (Z), with collision frequency: 



e 4 8-k 

V ee = v ee = n e — -2 In A (3.128) 

(47re ) 2 m\v\ 



v v ei = l -Zv ee (3.129) 



Hence (in the absence of other forces) 

d 



j (m e v) = - ( p u ee + p h> ei ) m p v (3.130) 

Z\ 

1 H ) v ee m e v (3.131) 



2 

This is equivalent to saying that the electron experiences an effective 'Frictional' force 

d ( Z \ 

F f = -n( m ev) = - [1 + — ) v ee m e v (3.132) 



dV V 2 

Z\ e 4 8tt In A 
2 J e (47re ) 2 m e v 2 



F f = 0)^7737^—^ (3.133) 
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Notice 

1. for Z = 1 slowing down is | on electrons | ions 

2. Fy decreases with i> increasing. 

Suppose now there is an electric field, E. The electron experiences an accelerating Force. 
Total force 

d , N / Z\ e 4 87rlnA . 

F = - raj) = -e.E + F f = -e£ - 1 + - ) n e ^ 3.134 

dt y 1 1 \ 2 J (47re ) 2 m e v 2 V ; 

Two Cases (When E is accelerating) 

1. \eE\ < \Ff\: Electron Slows Down 

2. \eE\ > \F f \: Electron Speeds Up! 

Once the electron energy exceeds a certain value its velocity increases continuously and the 
friction force becomes less and less effective. The electron is then said to ahve become a 



runaway . 
Condition: 



1 9 ( Z\ e 4 8vrlnA 

-m e v 2 > 1 + — )n e * — 3.135 

2 V 2 (47re ) 2 2eE V ; 



3.5.2 Plasma Resistivity (DC) 

Consider a bulk distribution of electrons in an electric field. They tend to be accelerated by 
E and decelerated by collisions. 

In this case, considering the electrons as a whole, no loss of total electron momentum by 
e — e collisions. Hence the friction force we need is just that due to V ei . 

If the electrons have a mean drift velocity Vd{« v t he) then 

^-(m e v d ) = -eE - V ei m e v d (3.136) 
at 

Hence in steady state 

v d = ^- . (3.137) 

m e v ei 

n e e 2 E 

j = -n e ev d = — — (3.138) 
m e v ei 

Now generally, for a conducting medium we define the conductivity, a, or resistivity, r], by 

j = aE ; m = E = ( 3 - 139 ) 
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The current is then 



Therefore, for a plasma, 

a = 1 = ^- (3.140) 

V m e u ei 

Substitute the value of V ei and we get 

n t Z 2 e 2 mf 8tt In A 
77 ~ • 5- (3.141) 

ne (47re ) 2 3y/2nTi 
Ze 2 ml 8tt In A 



(for a single ion species). (3.142) 

(47re )" 3V2n Tj 



3 

1 2 

- e 



Notice 

1. Density cancels out because more electrons means (a) more carriers but (b) more 
collisions. 

2. Main dependence is 77 oc T~ 3 / 2 . High electron temperature implies low resistivity (high 
conductivity). 

3. This expression is only approximate because the current tends to be carried by the 
more energetic electrons, which have smaller v e i] thus if we had done a proper average 
over f(v e ) we expect a lower numerical value. Detailed calculations give 

r7 = 5.2 x 1(T 5 — lnA 3 Vtm (3.143) 

{TjeVf* 

for Z — 1 (vs. ~ 1CT 4 in our expression). This is 'Spitzer' resistivity. The detailed 
calculation value is roughly a factor of two smaller than our calculation, which is not 
a negligible correction! 



3.5.3 Diffusion 

For motion parallel to a magnetic field if we take a typical electron, with velocity v\\ ~ v t e it 
will travel a distance approximately 

4 = v te /u ei (3.144) 

before being pitch-angle scattered enough to have its velocity randomised. [This is an order- 
of-magnitude calculation so we ignore P ee ] £ is the mean free path. 

Roughly speaking, any electron does a random walk along the field with step size £ and step 
frequency u ei . Thus the diffusion coefficient of this process is 

D e \\~e 2 e V ei ~ V ^- . (3.145) 

Similarly for ions 

A|| ^ £\vu ^ V S (3.146) 
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Notice 

vu/vei^(—) h ^- (if T e ~Ti) (3.147) 

Hence £ e ~ £i 

Mean free paths for electrons and ions are ~ same. 
The diffusion coefficients are in the ratio 

D fm \ I 

— - ~ I — - I : Ions diffuse slower in parallel direction. (3.148) 
Diffusion Perpendicular to Mag. Field is different 





Figure 3.7: Cross-field diffusion by collisions causing a jump in the gyrocenter (GC) position. 

Roughly speaking, if electron direction is changed by ~ 90° the Guiding Centre moves by 
a distance ~ tl- Hence we may think of this as a random walk with step size ~ and 
frequency V ei . Hence 

D el _ ~ r\ e v ei ~ V ^U ei (3.149) 

Ion transport is similar but requires a discussion of the effects of like and unlike collisions. 

Particle transport occurs only via unlike collisions. To show this we consider in more detail 
the change in guiding center position at a collision. Recall mv = qv A B which leads to 

vj_ = — r L A B (perp. velocity only). (3.150) 
m 

This gives 

FL = qB* (3 5 } 

At a collision the particle position does not change (instantaneously) but the guiding center 
position (r ) does. 

r' + r' L = r + r L Ar = r' - r = -(r' L - r L ) (3.152) 



60 



Change in r L is due to the momentum change caused by the collision: 

r ' L ~ TL = W 2 A m(v ' ~ v) " W 2 A A(mv) (3 ' 153) 

So 

Ar = -— A A(mv). (3.154) 

The total momentum conservation means that A(mv) for the two particles colliding is equal 
and opposite. Hence, from our equation, for like particles, Ar is equal and opposite. The 
mean position of guiding centers of two colliding like particles (r 01 +r 02 )/2 does not change. 

No net cross field particle (guiding center) shift. 

Unlike collisions (between particles of different charge q) do produce net transport of particles 
of either type. And indeed may move r i and r 2 in same direction if they have opposite 
charge. 

v 2 

D i± ~ r 2 Li Pu ie ~ ^ Pu ie (3.155) 
Notice that r 2 u /r\ e ~ mi/m e ; Pv~ ie /v ei ~ *± 

So D i± /D e ± ~ 1 (for equal temperatures). Collisional diffusion rates of particles are same 
for ions and electrons. 

However energy transport is different because it can occur by like-like collisins. 
Thermal Diffusivity: 

Xe ~ r\ e (V ei + I7 ee ) ~ r 2 Le V ei (V ei ~ I7 ee ) (3.156) 

Xi ~ r Li (^ie + V ii) - r Li V H i V U » V ie) (3.157) 



.2 T7 ™ ^2 



X,/Xe ~ ^^i^^ (equalT) (3.158) 

Collisional Thermal transport by Ions is greater than by electrons [factor ~ (mj/m e )^ ]. 
3.5.4 Energy Equilibration 

If T e 7^ Tj then there is an exchange of enegy between electrons and ions tending to make 
T e = Tj. As we saw earlier 

k 2m e m e ± 

Vei = P V ei = Vei (3.159) 



m, m 



So applying this to averages. 



K v ei ~ v ei (~ i/ ie ) (3.160) 



mi 
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Thermal energy exchange occurs ~ m e /mj slower than momentum exchange. (Allows T e ^ 
Ti). So 

<U] -^(Te-TO (3.161) 



e 

dt dt 



From this one can obtain the heat exchange rate (per unit volume), H ei , say: 

= ~jt (l n - T <) = i Gh T <) (3162) 

Q J Q 

= - i n-(T e -r i ) = -n^ ci (T e -T i ) (3.163) 



Important point: 

Zz/ ee ~ 

in; \m, 

'Electrons and Ions equilibrate among themselves much faster than with each other' 



3.6 Some Orders of Magnitude 

1. In A is very slowly varying. Typically has value ~ 12 to 16 for laboratory plasmas. 

2. V ei « 6 x Kr n (Vm 3 ) / (T e /e\/)^ (In A = 15, Z = 1). 

e.g. = 2 x 10 5 s" 1 (when n = 10 20 m~ 3 and T e = IkeV.) For phenomena which happen 
much faster than this, i.e. r<l/y e i~ 5/is, collisions can be ignored. 
Examples: Electromagnetic Waves with high frequency. 

3. Resistivity. Because most of the energy of a current carrying plasma is in the B field 
not the K.E. of electrons. Resistive decay of current can be much slower than z/ e j. E.g. 
Coaxial Plasma: (Unit length) 

Inductance L — [i In ^ 
Resistance R = rj 1 / na 2 
L/R decay time 

fi 7ia 2 b n e e 2 2 b 

t r ~ In - ~ — ^/x kcl m - 

77 a rripUpj a 



n e e 2 a 2 1 cu 2 a 2 1 1 
m e e c 2 F ei c 2 F ei z/ e 



» — . (3.165) 



Comparison 1 keV temperature plasma has ~ same (conductivity/) resistivity as a slab 
of copper (~2 x 10~ 8 Qm). 

Ohmic Heating Because rj oc T e 3 / 2 , if we try to heat a plasma Ohmically, i.e. by simply 
passing a current through it, this works well at low temperatures but its effectiveness falls 
off rapidly at high temperature. 
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Result for most Fusion schemes it looks as if Ohmic heating does not quite yet get us to the 
required ignition temperature. We need auxilliary heating, e.g. Neutral Beams. (These slow 
down by collisions.) 
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Chapter 4 

Fluid Description of Plasma 



The single particle approach gets to be horribly complicated, as we have seen. 

Basically we need a more statistical approach because we can't follow each particle separately. 
If the details of the distribution function in velocity space are important we have to stay 
with the Boltzmann equation. It is a kind of particle conservation equation. 



4.1 Particle Conservation (In 3-d Space) 



Ay 




Ax 

Figure 4.1: Elementary volume for particle conservation 

Number of particles in box AxAyAz is the volume, AV^ = AxAyAz, times the density n. 
Rate of change of number is is equal to the number flowing across the boundary per unit 
time, the flux. (In absence of sources.) 

d 

— -—-[AxAyAz n] = Flow Out across boundary. (4.1) 

L/C 

Take particle velocity to be v(r) [no random velocity, only flow] and origin at the center of 
the box refer to flux density as n\ = J. 

Flow Out = [J* (0,0, Az/2) - J z (0,0,-Az/2)] Ax Ay + x + y . (4.2) 
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Expand as Taylor series 

J z (0,0, V ) = J z (0) + ^-Jz-V (4.3) 
oz 



So, 



d 

flow out ~ — (nv z )AzAxAy + x + y (4.4) 
oz 



= AV V . (nv). 

Hence Particle Conservation 



d 

—n = V.(nv) (4.5) 



Notice we have essential proved an elementary form of Gauss's theorem 



V.Acfr = / A.dS. (4.6) 

Jd-y 

The expression: 'Fluid Description' refers to any simplified plasma treatment which does 
not keep track of v-dependence of / detail. 

1. Fluid Descriptions are essentially 3-d (r). 

2. Deal with quantities averaged over velocity space (e.g. density, mean velocity, ...). 

3. Omit some important physical processes (but describe others). 

4. Provide tractable approaches to many problems. 

5. Will occupy most of the rest of my lectures. 

Fluid Equations can be derived mathematically by taking moment^] of the Boltzmann Equa- 
tion. 

th moment / d 3 v (4.7) 
1st moment / vd 3 v (4-8) 
2nd moment J wd 3 v (4-9) 

These lead, respectively, to (0) Particle (1) Momentum (2) Energy conservation equations. 
We shall adopt a more direct 'physical' approach. 
1 They are therefore sometimes called 'Moment Equations.' 
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4.2 Fluid Motion 



The motion of a fluid is described by a vector velocity field v(r) (which is the mean velocity 
of all the individual particles which make up the fluid at r). Also the particle density 
n(r) is required. We are here discussing the motion of fluid of a single type of particle of 
mass/charge, m/q so the charge and mass density are qn and mn respectively. 

The particle conservation equation we already know. It is also sometimes called the 'Conti- 
nuity Equation' 

|-n + V.(nv) = (4.10) 
at 

It is also possible to expand the V. to get: 

d 

— n + (v.V)n + nV.v = (4.11) 
at 

The significance, here, is that the first two terms are the 'connective derivative" of n 

D d d , 



so the continuity equation can be written 

D 



n = -nV.v (4.13) 
Dt 



4.2.1 Lagrangian &; Eulerian Viewpoints 

There are essentially 2 views. 

1. Lagrangian . Sit on a fluid element and move with it as fluid moves. 




Figure 4.2: Lagrangean Viewpoint 



2. Eulerian . Sit at a fixed point in space and watch fluid move through your volume 
element: "identity" of fluid in volume continually changing 

means rate of change at fixed point (Euler). 

^ = J; = + v.V means rate of change at moving point (Lagrange). 

v.V = + + ■ change due to motion. 
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Figure 4.3: Eulerian Viewpoint 



Our derivation of continuity was Eulerian. From the Lagrangian view 



D _ d AN AN d _ 1 dAV 

Di U ~d~tAV = ~AV* Jt ~ ~ U AV ~W 



(4.14) 



since total number of particles in volume element (AN) is constant (we are moving with 
them). (AV = AxAyAz.) 



d 

Now — AV 
dt 



dAx . A dAy A . dAz . . 
AyAz + —^-AzAx + —^-AyAx 



dt 
AV 



dt 



dt 



1 dAx 1 dAy 1 dAz 

+ H 

Ax dt Ay dt Ax dt 



(4.15) 
(4.16) 



Hence 



and so 



But 



d(Ax) 
dt 



v x (Ax/2) - v x (-Ax/2) 

...y ...z 



A 9Vx 

Ax —— etc. 

ox 



dt 



AV = AV 



dv x | dv y | dv z 
dx dy dz 

D 



AVV.v 



Dt 



n = — nV.v 



(4.17) 
(4.18) 

(4.19) 
(4.20) 



Lagrangian Continuity. Naturally, this is the same equation as Eulerian when one puts 



+ v.V. 



D_ _ d_ 
Dt dt 

The quantity — V.v is the rate of (Volume) compression of element. 



4.2.2 Momentum (Conservation) Equation 

Each of the particles is acted on by the Lorentz force g[E + Uj AB] (uj is individual particle's 
velocity). 

Hence total force on the fluid element due to E-M fields is 



Y, (q [E + Ui A B]) = AN q (E + v A B) 



(4.21) 
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(Using mean: v = ^ u/AN.) 

E-M Force density (per unit volume) is: 

Fem = nq(E + v A B). (4.22) 

The total momentum of the element is 

^muj = m AN v = AV mnv (4.23) 

i 

so Momentum Density is mnv. 

If no other forces are acting then clearly the equation of motion requires us to set the time 
derivative of mnv equal to Fem- Because we want to retain the identity of the particles 
under consideration we want D/Dt i.e. the convective derivative (Lagrangian picture). 

In general there are additional forces acting. 

(1) Pressure (2) Collisional Friction. 



4.2.3 Pressure Force 

In a gas p(= nT) is the force per unit area arising from thermal motions. The surrounding 
fluid exerts this force on the element: 



Figure 4.4: Pressure forces on opposite faces of element. 



Net force in x direction is 

- p (Ax/2) AyAz + p (-Ax/2) AyAz (4.24) 

~ -AxAyAz ^ = -AV ^ = —AV (Vp) T (4.25) 
ox ox 

So (isotropic) pressure force density (/unit vol) 

F p = -Vp (4.26) 

How does this arise in our picture above? 

Answer : Exchange of momentum by particle thermal motion across the element boundary. 

Although in Lagrangian picture we move with the element (as defined by mean velocity v) 
individual particles also have thermal velocity so that the additional velocity they have is 

Wj = Uj — v 'peculiar' velocity (4.27) 
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Because of this, some cross the element boundary and exchange momentum with outside. 
(Even though there is no net change of number of particles in element.) Rate of exchange 

of momentum due to particles with peculiar velocity w, d 3 w across a surface element ds is 

/(w)raw d 3 w x yv_^ds (4.28) 
mom m density at w flow rate across ds 
Integrate over distrib function to obtain the total momentum exchange rate: 

ds. J mww/(w)(i 3 w (4.29) 

The thing in the integral is a tensor. Write 

p = J mww/(w)(i 3 w (Pressure Tensor) (4.30) 

Then momentum exchange rate is 

p . ds (4.31) 
Actually, if /(w) is isotropic (e.g. Maxwellian) then 

p xy = J m w x w y /(w)d 3 w = etc. (4.32) 

and p xx = J mwl f(w)d 3 w = nT{= p yy = p zz = 'p') (4.33) 

So then the exchange rate is pds. (Scalar Pressure). 

Integrate ds over the whole AV then x component of mom™ exchange rate is 

'Ax\ . . /— Ax~~ 

and so 

Total momentum loss rate due to exchange across the boundary per unit volume is 

Vp (= -F p ) (4.35) 

In terms of the momentum equation, either we put Vp on the momentum derivative side or 
F p on force side. The result is the same. 

Ignoring Collisions, Momentum Equation is 

^ (mnAVv) = [F EM + F p ] AV (4.36) 
Recall that nAV = AN ; jy t (AN) = 0; so 

D\ 

L.H.S. = mnAV — . (4.37) 
at 



P 



-) AyAz - -p (-^j AyAz = AV {Vp) x (4.34) 



Thus, substituting for F's: 

Momentum Equation. 



Dv fdv _ \ 

mn— = mn \ — + v. Vv = qn (E + v A B) - Vp (4.38) 
Dt \ at ) 
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4.2.4 Momentum Equation: Eulerian Viewpoint 

Fixed element in space. Plasma flows through it. 

1. E.M. force on element (per unit vol.) 

F EM = ng(E + v A B) as before. (4.39) 

2. Momentum flux across boundary (per unit vol) 

= V. J m(v + w)(v + w) /(w) d 3 w (4.40) 

= V.{Jm(vv + y w + wv +ww) /(w) rf 3 w} (4.41) 
integrates to 

= V.{mnvv + p} (4.42) 

= mn(v.V)t; + mv[V.(nv)] + Vj) (4.43) 

(Take isotropic p.) 

3. Rate of change of momentum within element (per unit vol) 

d 

= — (mnv) (4.44) 

Hence, total momentum balance: 
d 

— (mnv) + mn(v.V)v + mv [V.(nv)] + Vp = F EM (4.45) 

C6 

Use the continuity equation: 

^ + V.(nv) = , (4.46) 

to cancel the third term and part of the 1st: 

dn <9v <9v 

— (mnv) + mv (V. (nv)) = mv{— + V. (nv)} + mn— = mn— (4.47) 

Then take Vp to RHS to get final form: 
Momentum Equation: 



mn 



<9v 

w + (v.V)v 



nq (E + v A B) - Vp . (4.48) 
As before, via Lagrangian formulation. (Collisions have been ignored.) 
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4.2.5 Effect of Collisions 



First notice that like particle collisions do not change the total momentum (which is averaged 
over all particles of that species). 

Collisions between unlike particles do exchange momentum between the species. Therefore 
once we realize that any quasi-neutral plasma consists of at least two different species (elec- 
trons and ions) and hence two different interpenetrating fluids we may need to account for 
another momentum loss (gain) term. 

The rate of momentum density loss by species 1 colliding with species 2 is: 

vi 2 nimi(vi - v 2 ) (4.49) 
Hence we can immediately generalize the momentum equation to 
9vi 



mini 



+ (vi.V) vi 



dt 

With similar equation for species 2. 



nigi (E + vi A B) - Vpi - v 12 nimi (vi - v 2 ) (4.50) 



4.3 The Key Question for Momentum Equation: 

What do we take for pi 

Basically p = nT is determined by energy balance, which will tell how T varies. We could 
write an energy equation in the same way as momentum. However, this would then contain 
a term for heat flux, which would be unknown. In general, the k th moment equation contains 
a term which is a (k + l) th moment. 

Continuity, th equation contains v determined by 
Momentum, I s * equation contains p determined by 
Energy, 2 nd equation contains Q determined by . . . 

In order to get a sensible result we have to truncate this hierarchy. Do this by some sort of 
assumption about the heat flux. This will lead to an 
Equation of State: 

pn~"' = const. (4-51) 

The value of 7 to be taken depends on the heat flux assumption and on the isotropy (or 
otherwise) of the energy distribution. 

Examples 

1. Isothermal: T = const.: 7 = 1. 

2. Adiabatic/Isotropic: 3 degrees of freedom 7 = |. 

3. Adiabatic/1 degree of freedom 7 = 3. 
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4. Adiabatic/2 degrees of freedom 7 = 2. 



In general, n(£/2)5T = — p(SV/V) (Adiabatic £ degrees) 



£ 5T -SV 5n 

So 

* = *! + «: = ( 1+ |)*!, ,4.53) 
p n T \ £ / n 

i.e. 

p n -( 1+ !) = const. (4.54) 

In a normal gas, which 'holds together' by collisions, energy is rapidly shared between 3 
space-degrees of freedom. Plasmas are often rather collisionless so compression in 1 dimension 
often stays confined to 1-degree of freedom. Sometimes heat transport is so rapid that the 
isothermal approach is valid. It depends on the exact situation; so let's leave 7 undefined 
for now. 



4.4 Summary of Two-Fluid Equations 

Species j 
Plasma Response 



1. Continuity: 



dt 

2. Momentum: 



(> "' j + V.fav,-) = (4.55) 



= njqj (E + Vj A B) — Vpj — Vj k rijmj (vj — v k ) (4.56) 



3. Energy/Equation of State: 

Pj-nj" 1 = const.. (4.57) 

(j = electrons, ions). 



Maxwell's Equations 



V.B = V.E = p/e (4.58) 
1 Q"p f)TK 

VAB = wJ + ? - VAE=- (459) 
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With 



p = q e n e + qiiii = e (-n e + Zn^) 
j = g e n e v e + qjjiiVi = e (-n e v e + Zn^i 
= —en e (v e — Vjr) (Quasineutral) 



(4.60) 
(4.61) 
(4.62) 



Accounting 



Unknowns 

n e ,rii 2 



v e , v i 

E, B 



6 
2 
6 

16 



Equations 
Continuity e, i 
Momentum e, i 
State e, i 
Maxwell 



2 
6 
2 

8 

18 



frni 2 of Maxwell (V. equs) are redundant because can be deduced from others: e.g. 

d 



V.(VAE) = 
and V. (V A B) = 



dt 



(V.B) 



1 d 



+ V.E 



(4.63) 
(4.64) 



i- L dt 

So 16 equs for 16 unknowns. 

Equations still very difficult and complicated mostly because it is Nonlinear 

In some cases can get a tractable problem by 'linearizing'. That means, take some known 
equilibrium solution and suppose the deviation (perturbation) from it is small so we can 
retain only the 1st linear terms and not the others. 



4.5 Two-Fluid Equilibrium: Diamagnetic Current 

S^b: £ + I, £ = 0. 

Straight B-field: B = Bz. 
Equilibrium: g = (E = -V0) 

Collisionless: v — > 0. 
Momentum Equation(s): 

"W v r^) v / = ™j*(E + Vj A B) - Vpj (4.65) 
Drop j suffix for now. Then take x, y components: 

d dp 
mn Vx d~x Vx = Uq ^ + Vy ' ~ dx 

mn v x -^-v y = nq(0 — v x B) (4.67) 

(JjJb 
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Eg |4.67| is satisfied by taking v x = 0. Then [466 



i.e. 



or, in vector form: 



nq{E x + v y B) - ^ = 0. (4.68) 
ax 



~ Ex , 1 d P (A RO\ 

E A B Vp B fin . 

B A nq B l 

EAB drift Diamagnetic Drift 



Notice: 



• In magnetic field (_L) fluid velocity is determined by component of momentum equation 
orthogonal to it (and to B). 

• Additional drift (diamagnetic) arises in standard FAB form from pressure force. 

• Diagmagnetic drift is opposite for opposite signs of charge (electrons vs. ions). 

Now restore species distinctions and consider electrons plus single ion species i. Quasineu- 
trality says r^g; = —n e q e . Hence adding solutions 

n e q e v e + UiqiVi = + n e q e )^ -V (p e + p<) A — (4.71) 



Hence current density: 

j = -V(p e +p,)A-^ (4.72) 

This is the diamagnetic current. The electric field, E, disappears because of quasineutrality. 
(General case Yjj1j n jVj = — V (J2pj)/\B / B 2 ). 

4.6 Reduction of Fluid Approach to the Single Fluid 
Equations 

So far we have been using fluid equations which apply to electrons and ions separately. 
These are called 'Two Fluid' equations because we always have to keep track of both fluids 
separately. 

A further simplification is possible and useful sometimes by combining the electron and ion 
equations together to obtain equations governing the plasma viewed as a 'Single Fluid'. 
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Recall 2-fluid equations: 

UTl ' 

Continuity (Cj) + V.^v,) = 0. (4.73) 

Momentum (Mj) rrijUj ^ + Vj.vj v, = n^- (E + v, A B) - Vpj + F jk (4.74) 
(where we just write Fj k = —Vj k njrrij (vj — v fe ) for short.) 

Now we rearrange these 4 equations (2x2 species) by adding and subtracting appropriately 
to get new equations governing the new variables: 



Mass Density p m 


= n e m e + riifrii 


(4.75) 


C of M Velocity V 


= (n e m e v e + niiJiiVi) / p m 


(4.76) 


Charge density p q 


= q e n e + qiUi 


(4.77) 


Electric Current Density j 


= q e n e v e + qtriiVi 


(4.78) 




= q e n e (v e — Vj) by quasi neutrality 


(4.79) 


Total Pressure p 


= Pe+Pi 


(4.80) 



1st equation: take m e x C e + mi x Cj 



(1) ^ + V. (p m V) = Mass Conservation (4.81) 



2nd take q e x C e + qj x Cj 



(2) + V.j = Charge Conservation (4.82) 

3rd take M e + Mj. This is a bit more difficult. RHS becomes: 

(E + v, A B) - V pj + F jk = p q E + j A B - V (p e + Pi ) (4.83) 
(we use the fact that F ei — F ie so no net friction). LHS is 

' v,-vj r, (4.84) 

The difficulty here is that the convective term is non-linear and so does not easily lend itself 
to reexpression in terms of the new variables. But note that since m e « rrii the contribution 
from electron momentum is usually much less than that from ions. So we ignore it in this 
equation. To the same degree of approximation V ~ v«: the CM velocity is the ion velocity. 
Thus for the LHS of this momentum equation we take 

E m ^« f| + v rv) v, ~ Pm [| + V.V) V (4.85) 



3 
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so: 



(3) 



'd_ 

dt 



+ V. V V = pqE + j A B — Vp 



Finally we take ^M e + ^-M, to get: 
d 



n iQi 



dt 



+ (v,-.V) 



v,- = E{^(E + v,-AB)-^V ft - 



rrij 



Qj_ 

777-1 



F jk } 



(4.86) 



(4.87) 



Again several difficulties arise which it is not very profitable to deal with rigorously. Observe 
that the LHS can be written (using quasineutrality + n e q e = 0) as p m §i (^~) provided 
we discard the term in (v.V)v. (Think of this as a linearization of this question.) [The 
(v.V)v convective term is a term which is not satisfactorily dealt with in this approach to 
the single fluid equations.] 

In the R.H.S. we use quasineutrality again to write 



y n jQg 

^ rrii 



v, = 



2 2 



n e m e nim, 

2 

l ' Hq 

e 



E = niq 



2 2 mini + m e n e 



E 



q e qi 



n e m e nimi 



m e mi 



p m E, (4.88) 



777,; 



q e qi rn e q e mi n^me 

{ v e H Vj} 

m e mi q { q e 

q e qi 



q e qi 



m e mi 



I rrii m e 

-{n e m e \ e + UiiJiiWi 1 

V Qi Qe 

{p m v- — + — j} 
V Qi qe J 



(g e n e v e + qiUi-Vi)} 



Also, remembering F ei = -v ei n e m i {^ e - V;) = -F ie , 



777, 



-v e i n e q e - n e qi 



777, 



(v e - Vi) 



So we get 



dt \p„ 



q e qi 



m e mi 



p m E + <^ p m V 



777j 777 e 

v qi q e 



j AB 



Vp e V/7; 

777 e TTli 



q e m e \_. 

1 f«J 

9i rrii / 



(4.89) 



(4.90) 



(4.91) 



Regroup after multiplying by JL 
E+VAB = 



<?e<JiPm 

777 e 777i 8 ( j 



g e gi <9t Vp m / 

— Vj9 e + — Vp, 
777 P 777,- 



1 / 777j 


m e j . 


— — + 




Pm V Qi 


ft / 


\ m e m, 


•(- 







j AB 



q e m e 



m e mi 



qi rrii J q e qiPm 



Vei} 



(4.92) 
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Notice that this is an equation relating the Electric field in the frame moving with the fluid 
(L.H.S.) to things depending on current j i.e. this is a generalized form of Ohm's Law. 

One essentially never deals with this full generalized Ohm's law. Make some approximations 
recognizing the physical significance of the various R.H.S. terms. 

Tn e Tfli o ( j \ 

— - — -— — arises from electron inertia. 
q e qi dt \p m ) 

it will be negligible for low enough frequency. 

— ( — - H ] j A B is called the Hall Term. 

pm \q% q e J 

and arises because current flow in a B-field tends to be diverted across the magnetic field. 
It is also often dropped but the justification for doing so is less obvious physically. 

— Vpi term <C — Vp e for comparable pressures, 

mi m e 

and the latter is ~ the Hall term; so ignore (^Vpj/mj. 

Last term in j has a coefficient, ignoring m e /mi c.f. 1 which is 

m e miV ei m e V ei 

r = — — = i] the resistivity. (4.93) 

q e qi {nirrii) q l e n e 

Hence dropping electron inertia, Hall term and pressure, the simplified Ohm's law becomes: 

E + V A B = rjj (4.94) 

Final equation needed: state: 

p e n~ le + pin' 11 = constant. 
Take quasi- neutrality =>- n e oc n» oc p m . Take 7 e = 7«, then 

PPm 7 = const. (4.95) 

4.6.1 Summary of Single Fluid Equations: M.H.D. 

Mass Conservation : ^ + V (p m V) = (4.96) 

dp 

Charge Conservation : — — + V.j = (4.97) 

at 

Momentum : pm {jft + V ' V ) V = ^ E + J A B ~ V ^ ( 4 - 98 ) 

Ohm's Law : E + V A B = r] j (4.99) 

Eq. of State : pp~ 7 = const. (4.100) 
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4.6.2 Heuristic Derivation/ Explanation 

Mass Charge: Obvious. 

Mom m Pm ^ + V.vjv= p q E + jjVB - Vp (4.101) 

V „ ' Electric Magnetic Force Pressure 



rate of change of b ° d y f ° rc <= on current 

total momentum density 



Ohm's Law 



The electric field 'seen' by a moving (conducting) fluid is E + V A B = Ey electric field in 
frame in which fluid is at rest. This is equal to 'resistive' electric field rjy. 

Ey = E + VAB = ??j (4.102) 

The p q E term is generally dropped because it is much smaller than the j A B term. To see 
this, take orders of magnitude: 

V.E = p q /e so p q ~Ee /L (4.103) 

VAB = po3 (+^^J so vE = j~B/p L (4.104) 



Therefore 



p q E 6q ( B \ 2 Lp L 2 /c 2 f light transit time \ 



jB L \p aL J B 2 (p aL 2 ) 2 ^resistive skin time J ^ ^ 

This is generally a very small number. For example, even for a small cold plasma, say T e — 1 
eV (er m 2 x 10 3 mho/m), L = 1 cm, this ratio is about 10 -8 . 

Conclusion: the p q E force is much smaller than the j A B force for essentially all practical 
cases. Ignore it. 

Normally, also, one uses MHD only for low frequency phenomena, so the Maxwell displace- 
ment current, d'E/c 2 dt can be ignored. 

Also we shall not need Poisson's equation because that is taken care of by quasi-neutrality. 



4.6.3 Maxwell's Equations for MHD Use 

V.B = ; VAE=— ; V A B = /ij . (4.106) 

The MHD equations find their major use in studying macroscopic magnetic confinement 
problems. In Fusion we want somehow to confine the plasma pressure away from the walls 
of the chamber, using the magnetic field. In studying such problems MHD is the major tool. 

On the other hand if we focus on a small section of the plasma as we do when studying 
short-wavelength waves, other techniques: 2-fluid or kinetic are needed. Also, plasma is 
approx. uniform. 
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'Macroscopic' Phenomena MHD 
'Microscopic' Phenomena 2-Fluid/Kinetic 



4.7 MHD Equilibria 

Study of how plasma can be 'held' by magnetic field. Equilibrium =>- V = = 0. So 
equations reduce. Mass and Faraday's law are ~ automatic. We are left with 

(Mom m ) -> 'Force Balance' = j A B — Vp (4.107) 

Ampere V A B = (4.108) 

Plus V.B = 0, V.j = 0. 

Notice that provided we don't ask questions about Ohm's law. E doesn't come into MHD 
equilibrium. 

These deceptively simple looking equations are the subject of much of Fusion research. The 
hard part is taking into account complicated geometries. 

We can do some useful calculations on simple geometries. 



4.7.1 <9-pinch 



6_- |?^vcL 




Figure 4.5: #-pinch configuration. 

So called because plasma currents flow in ^-direction. 

Use MHD Equations 

Take to be oo length, uniform in z-dir. 

By symmetry B has only z component. 

By symmetry j has only 9- comp. 

By symmetry Vp has only r comp. 
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So we only need 



i.e. 



Force (j A B) r - (Vp) r = (4.109) 

Ampere (V A B) e = (/ij), (4.110) 

i.e. j 6 B z - ^P = (4.1H) 

8 B z = n j e (4.112) 



dr 

\i Q dr dr 



Eliminate j : _ * = ( 4 . 113 ) 



B 2 

Solution -^+p = const. (4.115) 
2/i G 




Figure 4.6: Balance of kinetic and magnetic pressure 

(4.116) 

[Recall Single Particle Problem] 

Think of these as a pressure equation. Equilibrium says total pressure = const. 

B 2 

+ p = const. (4.117) 

v — kinetic pressure 
magnetic pressure 
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Ratio of kinetic to magnetic pressure is plasma l (3\ 

P=^£ (4.H8) 

z 

measures 'efficiency' of plasma confinement by B. Want large (3 for fusion but limited by 
instabilities, etc. 

4.7.2 Z-pinch 




Figure 4.7: Z-pinch configuration, 
so called because j flows in z-direction. Again take to be oo length and uniform. 

j=j*e* B = B e e e (4.119) 



Eliminate j: 



or 



Force (j A B) r - (Vp) r = -j z B e - ^ = 



Ampere (V A B) 2 - (/x j) 
Bo d 



ld_ 

r dr 

dp 



[rBo) - fiojz = 



fi Q r dr 



(rBo) ~ £ = 



fi r dr \2jj, 



+ 



P) =0 



(4.120) 
(4.121) 

(4.122) 
(4.123) 



Extra Term Magnctic+Kinctic pressure 

Extra term acts like a magnetic tension force. Arises because B-field lines are curved. 
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Can integrate equation 



h Bl dr 



J a Ho 



+ 



fio r 



3. 

2/i 



+ p(r) 



If we choose b to be edge (p(b) = 0) and set a = r we get 




Figure 4.8: Radii of integration limits. 



(4 



2pL 2pL Jr Ho r' 



(4 



Force balance in z-pinch is somewhat more complicated because of the tension force 
can't choose p(r) and j(r) independently; they have to be self consistent. 

Example j = const. 



1 d . Hojz 

-— {rB e ) = fxojz =>• B e = — - r 
r or 2 



Hence 



p(r) 



1 /'/'.../ N 1 
2fi Q V 2 



{b 2 -r 2 + f b 2r'dr'} 

J r 



(4 

(4 
(4 



P 




Figure 4.9: Parabolic Pressure Profile. 

Also note B e (b) = Mfi so 

V= Wob- {h (4 
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4.7.3 'Stabilized Z-pinch' 

Also called 'screw pinch', 9 — z pinch or sometimes loosely just 'z-pinch'. 
Z-pinch with some additional B z as well as B e 

(Force), 3e B z - 3z B e - |- = (4.130) 
d 

Ampere : — B z = Li j e (4.131) 

— (4-132) 

Eliminate j: 

M-*° ( rB.)-% = * (4,33) 
/x or n r Or Or 

or 

— + tt (f + *) =° < 4 - 134 > 

Ma | ^" y sion Mag (0+z)+Kinetic pressure 

4.8 Some General Properties of MHD Equilibria 

4.8.1 Pressure & Tension 

j A B — Vp = : V A B = (4.135) 
We can eliminate j in the general case to get 

— (V A B) A B = Vp. (4.136) 

[J>o 

Expand the vector triple product: 

Vp = — (B.V) B —VB 2 (4.137) 

put b = |f| so that VB = VBb = BVb + bVB. Then 

Vp = — {B 2 (b.V)b + Bb(b.V)B}-^-VB 2 (4.138) 

fX a ^fXo 

= — (b.V)b- — (V-b(b.V))5 2 (4.139) 
fx 2/i G 

B 2 I B 2 \ 

= -(b.V)b-Vj— (4.140) 
\± y zfXp j 
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Now V_i_ (|^J is the perpendicular (to B) derivative of magnetic pressure and (b.V)b is the 
curvature of the magnetic field line giving tension. 

| (b.V)b | has value ^. i?: radius of curvature. 

4.8.2 Magnetic Surfaces 

= B. [j A B — Vp] = -B.Vp (4.141) 

*Pressure is constant on a field line (in MHD situation). 
(Similarly, = j.[j A B — Vp] = j.Vp.) 




p= const, contours 



Figure 4.10: Contours of pressure. 

Consider some arbitrary volume in which Vp ^ 0. That is, some plasma of whatever shape. 
Draw contours (surfaces in 3-d) on which p = const. At any point on such an isoberic surface 
Vp is perp to the surface. But B.Vp = implies that B is also perp to Vp. 




Figure 4.11: B is perpendicular to Vp and so lies in the isobaric surface. 

Hence 

B lies in the surface p = const. 

In equilibrium isobaric surfaces are 'magnetic surfaces'. 

[This argument does not work if p = const, i.e. Vp = 0. Then there need be no magnetic 
surfaces.] 
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4.8.3 'Current Surfaces' 



Since j.Vp = in equilibrium the same argument applies to current density. That is 

j lies in the surface p = const. 

Isobaric Surfaces are 'Current Surfaces'. 

Moreover it is clear that 

'Magnetic Surfaces' are 'Current Surfaces'. 

(since both coincide with isobaric surfaces.) 

[It is important to note that the existence of magnetic surfaces is guaranteed only in the 
MHD approximation when Vp 7^ > Taking account of corrections to MHD we may not 
have magnetic surfaces even if Vp 7^ 0.] 



4.8.4 Low (3 equilibria: Force-Free Plasmas 

In many cases the ratio of kinetic to magnetic pressure is small, /3 « 1 and we can approx- 
imately ignore Vp. Such an equilibrium is called 'force free'. 

j A B = (4.142) 

implies j and B are parallel, 
i.e. 

j = /x(r)B (4.143) 
Current flows along field lines not across. Take divergence: 

= V.j = V. (ji (r) B) = n (r) V.B + (B. V) fi (4.144) 
= (B.V)/x. (4.145) 

The ratio j/B = fi is constant along field lines. 

H is constant on a magnetic surface. If there are no surfaces, fi is constant everywhere. 



Example: Force-Free Cylindrical Equil. 



j a B = j = (i(r)B (4.146) 

VAB = fi j = fi of i{r)B (4.147) 

This is a somewhat more convenient form because it is linear in B (for specified fJ>(r)). 

Constant-^ : V A B = /i /iB (4.148) 

leads to a Bessel function solution 
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B z = B J (fi fjir) 
B e = BoJ^/io/ir) 



(4.149) 
(4.150) 



for /2 jj,r > 1st zero of J Q the toroidal field reverses. There are plasma confinement schemes 
with /i ~ const. 'Reversed Field Pinch'. 

4.9 Toroidal Equilibrium 

Bend a z-pinch into a torus 




Figure 4.12: Toroidal z-pinch 

Bq fields due to current are stronger at small R side =>- Pressure (Magnetic) Force outwards. 
Have to balance this by applying a vertical field B„ to push plasma back by A T5 V . 




Figure 4.13: The field of a toroidal loop is not an MHD equilibrium. Need to add a vertical 
field. 

Bend a #-pinch into a torus: B<p is stronger at small R side =^ outward force. 

Cannot be balanced by B„ because no j<f,. No equilibrium for a toroidally symmetric #-pinch. 

Underlying Single Particle reason: 

Toroidal #-pinch has B^ only. As we have seen before, curvature drifts are uncompensated 
in such a configuration and lead to rapid outward motion. 
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Figure 4.14: Charge-separation giving outward drift is equivalent to the lack of MHD toroidal 
force balance. 



We know how to solve this: Rotational Transform: get some B e . Easiest way: add j^. From 
MHD viewpoint this allows you to push the plasma back by A force. Essentially, this 
is Tokamak. 



4.10 Plasma Dynamics (MHD) 



When we want to analyze non-equilibrium situations we must retain the momentum terms. 
This will give a dynamic problem. Before doing this, though, let us analyse some purely 
Kinematic Effects. 

'Ideal MHD 7 Set eta = in Ohm's Law. 

A good approximation for high frequencies, i.e. times shorter than resistive decay time. 

E + VAB = 0. Ideal Ohm's Law. (4.151) 

Also 

V A E = Faraday's Law. (4.152) 

Together these two equations imply constraints on how the magnetic field can change with 
time: Eliminate E: 

rfYK 

+ VA(VAB) = +— (4.153) 
This shows that the changes in B are completely determined by the flow, V. 



4.11 Flux Conservation 

Consider an arbitrary closed contour C and spawning surface S in the fluid. 
Flux linked by C is 

$ = f B.ds (4.154) 
J s 

Let C and S move with fluid: 

Total rate of change of <£> is given by two terms: 
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Contour at time t, t'=t+dt 

Figure 4.15: Motion of contour with fluid gives convective flux derivative term. 



OB 



.ds 



+ f B.(VAdl) 



c 



Due to changes in B Due to motion of C 

- / VAE.ds- /(VAB).dl 

Js Jc 

— I (E + V A B).dl = by Ideal Ohm's Law. 

J c 



Flux through any surface moving with fluid is conserved. 



(4.155) 

(4.156) 
(4.157) 



4.12 Field Line Motion 

Think of a field line as the intersection of two surfaces both tangential to the field everywhere: 




Field Line 



Figure 4.16: Field line defined by intersection of two flux surfaces tangential to field. 
Let surfaces move with fluid. 

Since all parts of surfaces had zero flux crossing at start, they also have zero after, (by flux 
conserv.). 



Surfaces are tangent after motion 



=>- Their intersection defines a field line after. 

We think of the new field line as the same line as the old one (only moved). 
Thus: 

1. Number of field lines (= flux) through any surface is constant. (Flux Cons.) 

2. A line of fluid that starts as a field line remains one. 



4.13 MHD Stability 

The fact that one can find an MHD equilibrium (e.g. z-pinch) does not guarantee a useful 
confinement scheme because the equil. might be unstable. Ball on hill analogies: 




Stable Unstable Unstable 



Figure 4.17: Potential energy curves 

An equilibrium is unstable if the curvature of the 'Potential energy surface' is downward 
away from equil. That is if ^2{W pot } < 0. 

In MHD the potential energy is Magnetic + Kinetic Pressure (usually mostly magnetic). 

If we can find any type of perturbation which lowers the potential energy then the equil is 
unstable. It will not remain but will rapidly be lost. 

Example Z-pinch 

We know that there is an equilibrium: Is it stable? 
Consider a perturbation thus: 











1 1 

J \ 

1 




J 






■ i 


s 


■ ^ J 


















I j 











Figure 4.18: 'Sausage' instability 



Simplify the picture by taking the current all to flow in a skin. We know that the pressure 
is supported by the combination of B 2 /2fj, Q pressure and ~ tension forces. 



89 





j 




_p_ 








>— 


, B 



r 



Figure 4.19: Skin-current, sharp boundary pinch. 
At the place where it pinches in (A) 

Be and - increase — > Mag. pressure & tension increase =>• inward force no longer balance 
by p =>- perturbation grows. 

At place where it bulges out (B) 

Bg & - decrease Pressure & tension perturbation grows. 

Conclusion a small perturbation induces a force tending to increase itself. Unstable (= < 
0). 

4.14 General Perturbations of Cylindrical Equil. 

Look for things which go like exp[i(kz + md)\. [Fourier (Normal Mode) Analysis]. 




Figure 4.20: Types of kink perturbation. 
Generally Helical in form (like a screw thread). Example: m = 1 k ^ z-pinch 

4.15 General Principles Governing Instabilities 

(1) They try not to bend field lines. (Because bending takes energy). Perturbation 
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Figure 4.21: Driving force of a kink. Net force tends to increase perturbation. Unstable. 




^3 . 



Figure 4.22: Alignment of perturbation and field line minimizes bending energy. 

(Constant surfaces) lies along magnetic field. 
Example: #-pinch type plasma column: 




Figure 4.23: 'Flute' or 'Interchange' modes. 

Preferred Perturbations are 'Flutes' as per Greek columns — > 'Flute Instability' [Better 
name: 'Interchange Instability', arises from idea that plasma and vacuum change places.] 

(2) Occur when a 'heavier' fluid is supported by a 'lighter' (Gravitational analogy). 

Why does water fall out of an inverted glass? Air pressure could sustain it but does not 
because of Ray leigh- Taylor instability. 

Similar for supporting a plasma by mag field. 

(3) Occur when | B | decreases away from the plasma region. 
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© m %-^a\A 

Figure 4.24: Inverted water glass analogy. Rayleigh Taylor instability. 



R2 R2 

2ji a 2fi Q 



(4.158) 



Perturbation Grows. 



(4) Occur when field line curvature is towards the plasma (Equivalent to (3) because of 
V A B = in a vacuum) . 




O O £> <P 

Figure 4.25: Vertical upward field gradient is unstable. 
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5 ' 'E^U tt,n*r 





Figure 4.26: Examples of magnetic configurations with good and bad curvature. 
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4.16 Quick and Simple Analysis of Pinches 



#-pinch | B | = const, outside pinch 

= No field line curvature. Neutral stability 

z-pinch V | B | away from plasma outside 

= Bad Curvature (Towards plasma) =>- Instability. 

Generally it is difficult to get the curvature to be good everywhere. Often it is sufficient 
to make it good on average on a field line. This is referred to as 'Average Minimum B'. 
Tokamak has this. 

General idea is that if field line is only in bad curvature over part of its length then to perturb 
in that region and not in the good region requires field line bending: 



^ Perturbation prefered by 
--curvature, bends field line. 



Plasma Region 



Good Curvature 



Figure 4.27: Parallel localization of perturbation requires bending. 

But bending is not preferred. So this may stabilize. 

Possible way to stabilize configuration with bad curvature: Shear 
Shear of Field Lines 




Figure 4.28: Depiction of field shear. 



Direction of B changes. A perturbation along B at z 3 is not along B at z 2 or z\ so it would 
have to bend field there — > Stabilizing effect. 

General Principle: Field line bending is stabilizing. 

Example: Stabilized z-pinch 

Perturbations (e.g. sausage or kink) bend B z so the tension in B z acts as a restoring force 
to prevent instability. If wave length very long bending is less. =>- Least stable tends to be 
longest wave length. 

Example: 'Cylindrical Tokamak' 
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Tokamak is in some ways like a periodic cylindrical stabilized pinch. Longest allowable wave 
length = 1 turn round torus the long way, i.e. 



kR=l: 



X = ItxR. 



(4.159) 



Express this in terms of a toroidal mode number, n (s.t. perturbation oc expi(n0 + mO): 
<p = -| n — kR. 

Most unstable mode tends to be n — 1. 

[Careful! Tokamak has important toroidal effects and some modes can be localized in the 
bad curvature region (n^ 1). 



Figure 4.29: Ballooning modes are localized in the outboard, bad curvature region. 



/ 



Good 
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Chapter 5 



Electromagnetic Waves in Plasmas 



5.1 General Treatment of Linear Waves in Anisotropic 
Medium 

Start with general approach to waves in a linear Medium: Maxwell: 

1 <9E <9B 

we keep all the medium's response explicit in j. Plasma is (infinite and) uniform so we Fourier 
analyze in space and time. That is we seek a solution in which all variables go like 

expi(k.x — uot) [real part of] (5.2) 

It is really the linearised equations which we treat this way; if there is some equilibrium field 
OK but the equations above mean implicitly the perturbations B, E, j, etc. 

Fourier analyzed: 

ik A B = /ij + — ^E ; ik A E = iuB (5.3) 

cr 

Eliminate B by taking kA second eq. and wx 1st 

ioo 2 

ik A (k A E) = LUfioj -E (5.4) 

cr 



So 

k A (k A E) + ^-E + iufi i = (5.5) 



cu 2 



c 2 



Now, in order to get further we must have some relationship between j and E(k, u). This 
will have to come from solving the plasma equations but for now we can just write the most 
general linear relationship j and E as 

j = o-.E (5.6) 
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cr is the 'conductivity tensor'. Think of this equation as a matrix e.g.: 

(5.7) 





This is a general form of Ohm's Law. Of course if the plasma (medium) is isotropic (same 
in all directions) all off-diagonal a's are zero and one gets j = cxE. 

Thus 



2 



k(k.E) — k 2 E + ^-E + iujfi cr.E = (5.8) 
cr 

Recall that in elementary E&M, dielectric media are discussed in terms of a dielectric con- 
stant e and a "polarization" of the medium, P, caused by modification of atoms. Then 

e D E = E^ - ^P, and V.D = p cxt (5.9) 

Displacement Polarization externalcharge 

and one writes 

P= ^ e E (5.10) 
susceptibility 

Our case is completely analogous, except we have chosen to express the response of the 
medium in terms of current density, j, rather than "polarization" P For such a dielectric 
medium, Ampere's law would be written: 

1 <9D d 

V A B = j ext + — = -ee E, if j ext = , (5.11) 
fx at at 

where the dielectric constant would be e = 1 + x- 

Thus, the explicit polarization current can be expressed in the form of an equivalent dielectric 
expression if 

dE „ <9E d . n , 

J + e - = *.E + e- = -e e.E (5.12) 

or 

e = l + — — (5.13) 
— iue 

Notice the dielectric constant is a tensor because of anisotropy. The last two terms come 
from the RHS of Ampere's law: 

j + |(e E) . (5.14) 

If we were thinking in terms of a dielectric medium with no explicit currents, only implicit (in 
e) we would write this ^(ee D E); e the dielectric constant. Our medium is possibly anisotropic 
so we need J^(e D e.E) dielectric tensor. The obvious thing is therefore to define 

e = 1 + — L_ a- = i + (5.15) 

— lUJto UJ 
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Then 



id 



k(k.E) - k 2 E + — e.E = 



and we may regard e(k, 00) as the dielectric tensor. 
Write the equation as a tensor multiplying E: 



D.E = 



with 



D = {kk-A; 2 l + — e} 



Again this is a matrix equation i.e. 3 simultaneous homogeneous eqs. for E. 



D X X D X y 



D 



yx 



D, 



( E x 

Ey 

\ E z 



In order to have a non-zero E solution we must have 

det I D 1= 0. 



(5.16) 



(5.17) 
(5.18) 

(5.19) 



(5.20) 



This will give us an equation relating k and u, which tells us about the possible wavelengths 
and frequencies of waves in our plasma. 

5.1.1 Simple Case. Isotropic Medium 



<T = a 1 
e = e 1 

Take k in z direction then write out the Dispersion tensor D. 



D 















: .( 


V 





kk / I 



kk 



Poo 

k 2 
k 2 

k 2 l 



+ 



4e 



\ 

V C ) 





J 1 
-pi' 



(5.21) 
(5.22) 



-k 2 + %e 





Take determinant: 

det I D |= ^-h 
Two possible types of solution to this dispersion relation: 











k 2 + 


c z 










4e 


cu 2 


\ 2 u 2 


e = 




j ^ 



(5.23) 



(5.24) 
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(A) 



k 2 + —e 



= 0. 





















s ) 



















E z = 0. 



(5.25) 
(5.26) 



Electric field is transverse (E.k = 0) 
Phase velocity of the wave is 



k 



c 

7^ 



This is just like a regular EM wave traveling in a medium with refractive index 

kc 



N 



(B) 





^ D xx 




















V o 









-e = i.e. e = 



Er.. 



Ey \ = 

E z 



E x — E y — 0. 



(5.27) 

(5.28) 

(5.29) 
(5.30) 



Electric Field is Longitudinal (E A k = 0) E || k. 

This has no obvious counterpart in optics etc. because e is not usually zero. In plasmas 
e = is a relevant solution. Plasmas can support longitudinal waves. 



5.1.2 General Case (k in z-direction) 



D 



-N + e xx 

£y X 



x y 

-N 2 + e 

e zy 



N 2 



k 2 c 



(5.31) 



When we take determinant we shall get a quadratic in iV 2 (for given uS) provided e is not 
explicitly dependent on k. So for any uo there are two values of iV 2 . Two 'modes'. The 
polarization E of these modes will be in general partly longitudinal and partly transverse. 
The point: separation into distinct longitudinal and transverse modes is not possible in 
anisotropic media (e.g. plasma with B Q ). 

All we have said applies to general linear medium (crystal, glass, dielectric, plasma). Now 
we have to get the correct expression for cr and hence e by analysis of the plasma (fluid) 
equations. 
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5.2 High Frequency Plasma Conductivity 



We want, now, to calculate the current for given (Fourier) electric field E(k, u), to get the 
conductivity, cr. It won't be the same as the DC conductivity which we calculated before 
(for collisions) because the inertia of the species will be important. In fact, provided 

to > 9 ei (5.32) 

we can ignore collisions altogether. Do this for simplicity, although this approach can be 
generalized. 

Also, under many circumstances we can ignore the pressure force — Vp. In general will 
be true if | 3> v te ,i We take the plasma equilibrium to be at rest: v Q = 0. This gives a 
manageable problem with wide applicability. 

Approximations: 

Collisionless u ei = 

'Cold Plasma' Vp = (e.g. T ~ 0) (5.33) 

Stationary Equil v D = 

5.2.1 Zero B-field case 

To start with take B D = 0: Plasma isotropic Momentum equation (for electrons first) 

<9v 



ran 



at + (v ' v)v 



nqE (5.34) 



Notice the characteristic of the cold plasma approx. that we can cancel n from this equation 
and on linearizing get essentially the single particle equation. 

m ~0^ = <?E (Drop the 1 suffix now). (5.35) 
This can be solved for given to as 

v = — - — E (5.36) 
—turn 

and the current (due to this species, electrons) is 

2 

j = ngv = — — ^ — E (5.37) 
— lum 

So the conductivity is 

(T = i 7 ^- (5.38) 
uim 

Hence dielectric constant is 

e = l + — * = !-(— )- 2 =l + X (5.39) 
ooe \me ) oo z 
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Longitudinal Waves (B G = 0) 

Dispersion relation we know is 

e = »=l-Wi (5.40) 
\me J uj z 

[Strictly, the e we want here is the total e including both electron and ion contributions to 
the conductivity. But 

^ ~ ^ (for z = 1) (5.41) 

so to a first approximation, ignore ion motions.] 
Solution 

^=(^-). (5.42) 
In this approx. longitudinal oscillations of the electron fluid have a single unique frequency: 

^69'- (5 - 43) 

This is called the 'Plasma Frequency ' (more properly u pe the 'electron' plasma frequency). 
If we allow for ion motions we get an ion conductivity 



iniqf 
o~i = 

(jJUli 

and hence 



(5.44) 



e to t = H \<Je + Oi) = 1 1 — (5.45) 



= 1 - (^L + / 

where 



uj 2 



is the 'Ion Plasma Frequency'. 



Simple Derivation of Plasma Oscillations 

Take ions stationary; perturb a slab of plasma by shifting electrons a distance x. Charge 
built up is n e qx per unit area. Hence electric field generated 

E = - 7 ^- (5.47) 
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Figure 5.1: Slab derivation of plasma oscillations 



Equation of motion of electrons 

dv n e q 2 e x 

m e — = ; (5.48) 

at t Q 

Simple harmonic oscillator with frequency 

(n e q 2 \ ^ 
e e Plasma Frequency. (5.50) 

e m e J 

The Characteristic Frequency of Longitudinal Oscillations in a plasma. Notice 

1. uj = uj p for all k in this approx. 

2. Phase velocity | can have any value. 

3. Group velocity of wave, which is the velocity at which information/energy travel is 

».= £ = (.,. (5.51) 

In a way, these oscillations can hardly be thought of as a 'proper' wave because they do 
not transport energy or information. (In Cold Plasma Limit). [Nevertheless they do emerge 
from the wave analysis and with less restrictive approxs do have finite v g ] 

Transverse Waves (B Q = 0) 
Dispersion relation: 

uj 2 

- k 2 + — e = (5.52) 
cr 
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or 



N 2 EE 



k 2 c 2 



i - (u 2 pe + /u 2 



CO 



/co 2 



(5.53) 




Figure 5.2: Unmagnetized plasma transverse wave. 




Figure 5.3: Alternative dispersion plot. 



Alternative expression: 



-* 2 + ^-4 = (5.54) 
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which implies 



2 2 , ; 2 2 

OJ = OJZ + fee 



= (*;+*v) 



(5.55) 
(5.56) 



5.2.2 Meaning of Negative N 2 : Cut Off 

When A^ 2 < (for uo < uo p ) this means iV is pure imaginary and hence so is k for real uo. 
Thus the wave we have found goes like 



exp{± | k | x — iut} 



(5.57) 



i.e. its space dependence is exponential not oscillatory. Such a wave is said to be 'Evanescent' 
or 'Cut Off'. It does not truly propagate through the medium but just damps exponentially. 

Example: 




Figure 5.4: Wave behaviour at cut-off. 



A wave incident on a plasma with > u 2 is simply reflected, no energy is transmitted 
through the plasma. 



5.3 Cold Plasma Waves (Magnetized Plasma) 

Objective: calculate e,D,k(u;), using known plasma equations. 
Approximation: Ignore thermal motion of particles. 

Applicability: Most situations where (1) plasma pressure and (2) absorption are negligible. 
Generally requires wave phase velocity ^> ^thermal- 
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5.3.1 Derivation of Dispersion Relation 

Can "derive" the cold plasma approx from fluid plasma equations. Simpler just to say that 
all particles (of a specific species) just move together obeying Newton's 2nd law: 



<9v 

m— = qCE + v A B) 
ot 



(5.58) 



Take the background plasma to have E = 0, B = B and zero velocity. Then all motion 
is due to the wave and also the wave's magnetic field can be ignored provided the particle 
speed stays small. (This is a linearization). 



<9v 

m— =g(E + vAB ), 



(5.59) 



where v, E oc expi(k.x — out) are wave quantities. 

at 



Substitute — > — iuo and write out equations. Choose axes such that B = -Bq(0,0, 1). 



— luumvx 



q{E x + v y B ) 
q(E y - v x B Q ) 



(5.60) 



Solve for v in terms of E. 



iumv z = qE z 

q ( iuE x — QE y 



v x = 



m 



00" 



n 2 



v, = 



q / VtE x + iuEy 

m\ uj 2 — Vt 2 
q i 



(5.61) 



-E z 



m uo 



where Vt = ^ is the gyrofrequency but its sign is that of the charge on the particle species 
under consideration. 

Since the current is j = qvn = cr.E we can identify the conductivity tensor for the species 
(j) as: 





iw 








rrij 




rrij u> 2 —Q'j 


rrij 




rrij u) 2 —Q 2 
J j 










m u> 



(5.62) 



The total conductivity, due to all species, is the sum of the conductivities for each 



(5.63) 



So 



'yy 



q{rij luj 
f rrij uo 2 — Vt 2 

3 J 



E 



(5.64) 
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a 



xy 



—a 



yx 



-E 



Tj n 3 



rrij uj 2 — Q 2 
3 j 



(7, 



<lj n 3 i 

j m 3 U 



Susceptibility x — z^— a - 



^xx ^xy 




' S 


-iD 





£yx £yy 




iD 


S 





e zz _ 










p 



where 



-yy 



-ze 



5=1-5: 



CO 



U P3 

2 -n 2 



yx 



= jj = ^2^j 



UJ UJ C 



UJ 



UJ 



and 



UJ 



Qj n 3 

e rrij 



is the "plasma frequency" for that species. 
S & D stand for "Sum" and "Difference" : 

S=\{R + L) D= l -{R-L) 

where R & L stand for "Right-hand" and "Left-hand" and are: 



UJ 



P'J 



uj (uj + fL 



i-E 



UJ. 



p.) 



uj (uj — Qj 



(5.65) 
(5.66) 



(5.67) 



(5.68) 
(5.69) 
(5.70) 

(5.71) 



(5.72) 



(5.73) 



The R & L terms arise in a derivation based on expressing the field in terms of rotating 
polarizations (right & left) rather than the direct Cartesian approach. 

We now have the dielectric tensor from which to obtain the dispersion relation and solve 
it to get k(u;) and the polarization. Notice, first, that e is indeed independent of k so the 
dispersion relation (for given a;) is a quadratic in N 2 (or k 2 ). 

Choose convenient axes such that k y = N y = 0. Let 9 be angle between k and B so that 

N z = Ncos9, N x = Nsin9. (5.74) 

Then 



D = 



-N 2 cos 2 9 + S -iD N 2 sin 9 cos 9 
+iD -N 2 + S 

N 2 sin 9 cos 9 -N 2 sin 2 9 + P 



(5.75) 
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and 
where 



Solutions are 



|| D ||= AN*-BN 2 + C (5.76) 

A = Ssm 2 6 + Pcos 2 6 (5.77) 

B = RLsm 2 9 + PS(l + cos 2 #) (5.78) 

C = PRL (5.79) 



N 2 = ^f, (5.80) 



where the discriminant, F, is given by 

F 2 = (RL - PS) 2 sin 4 9 + AP 2 D 2 cos 2 6 (5.81) 

after some algebra. This is often, for historical reasons, written in the equivalent form (called 
the Appleton-Hartree dispersion relation) 

, 2(A-B + C) , 

N 2 = 1 - 4i J 5.82 

2A - B ± F y ' 

The quantity F 2 is generally +ve, so N 2 is real =>- "propagating" or "evanescent" no wave 
absorption for cold plasma. 

Solution can also be written 

2 „ P (N 2 - R) (N 2 - L) 

This compact form makes it easy to identify the dispersion relation at 6 = & | i.e. parallel 
and perpendicular propagation tan# = 0, oo. 

Parallel: P = , N 2 = R N 2 = L 

Perp: N 2 = ^ N 2 = P . 

Example: Right-hand wave 

iV 2 = R. (Single Ion Species). 

2 2 

N 2 = l~ , Upe lr , „ - , „ (5.84) 

u (u - \Q, e \) uj(uj + \Vti\) v ' 

This has a wave resonance N 2 — > oo at a; = |il e |, only. Right-hand wave also has a cutoff at 
i? = 0, whose solution proves to be 



| ^\ g | | | 

UJ = LU R = h 



i \ 2 



2 I + ^pe + ^pi 



1/2 

(5.85) 
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Since to, ^> m e this can be approximated as: 

i a 



Mr ^ 



1+1 + 4 



pe 



This is always above \£l e 




(5.86) 



Figure 5.5: The form of the dispersion relation for RH wave. 

One can similarly investigate LH wave and perp propagating waves. The resulting wave 
resonances and cut-offs depend only upon 2 properties (for specified ion mass) (1) Density 
<-> u 2 e (2) Magnetic Field <-> \Q e \. [Ion values u P i, are got by ^ factors.] 

These resonances and cutoffs are often plotted on a 2-D plane ^ , ^| (oc B, n) called the 

C M A Diagram. 

We don't have time for it here. 



5.3.2 Hybrid Resonances Perpendicular Propagation 

"Extraordinary" wave N 2 = ^ 



N 2 



{u + a) (u + a) - ^ (u + a) - ^ (w + a 



(w-fie) (oj - Qi) - ^ (u - Qi) ... 



Resonance is where denominator = 0. Solve the quadratic in u 2 and one gets 



u 2 = 



o& + n| + a£+jg ± 



a;, 



pe 



fi 2 



a 2 



(5.87) 



(5.88) 
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Neglecting terms of order ^ (e.g. U ^ L ) one gets solutions 



u UH = uj pe + Q e Upper Hybrid Resonance. 



LH 



e pi 



9? p + uj 



2 

pe 



Lower Hybrid Resonance. 



At very high density, u 2 e 2> fl 



u 2 LH ~ \n e \\tti\ 



geometric mean of cyclotron frequencies. 
At very low density, uj 2 e <C Vt 2 e 

ion plasma frequency 

Usually in tokamaks u) 2 e ~ f2 2 . Intermediate. 
Summary Graph (fi > uo p ) 



(5.89) 
(5.90) 

(5.91) 
(5.92) 




Figure 5.6: Summary of magnetized dispersion relation 

Cut-offs are where A^ 2 = 0. 
Resonances are where A^ 2 — > oo. 

Intermediate angles of propagation have refractive indices between the # = 0, | lines, in the 
shaded areas. 
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5.3.3 Whistlers 



(Ref. R.A. Helliwell, "Whistlers & Related Ionospheric Phenomena," Stanford UP 1965.) 
For N 2 ^> 1 the right hand wave can be written 

,2 



Group velocity is 



Then since 



we have 



duo 



(Nuo) 



Thus 



Group Delay is 



N 2 



—oj. 



pe 



UO {uj — |Q e |) 

duo I dk\ 1 
' Jg = dk = \~Io~ ) 

N = 



(N = kc/uo) 



d /Nu\ 
duo V c / 



n -l 



UOn 



uo 2 (|Q e | — a;) 



da; ( 


Q e \ - uo)' 2 




Up/ 2 


(|fl e 


\ 5 1 

— UO) 2 U02 




V 1 ^ e 1 




\^ 1 

— C<j) 2 U02 




c 



UOr, 



1 ) 

2 



c 2(|n c | -uo) 2 uo2 



uo p \Vt e \ 



;U0 2 



— oc — 3 oc 1 

V 9 002 ( Q e — tS) 2 (j!L.Y(-\-JiL.\ 



(5.93) 

(5.94) 
(5.95) 



(5.96) 

(5.97) 
(5.98) 




Figure 5.7: Whistler delay plot 
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Plot with — as x-axis. 
Resulting form explains downward whistle. 
Lightning strike ~ 5-function excites all frequencies. 
Lower ones arrive later. 



Examples of actual whistler sounds can be obtained from http://www-istp.gsfc.nasa 



gov/ istp/polar/polar_pwi_sounds .html 



5.4 Thermal Effects on Plasma Waves 

The cold plasma approx is only good for high frequency, N 2 ~ 1 waves. If u is low or N 2 ^> 1 
one may have to consider thermal effects. From the fluid viewpoint, this means pressure. 
Write down the momentum equation. (We shall go back to Bq = 0) linearized 

mn—- = nqEii — Vpi ; (5.99) 
remember these are the perturbations: 

P = Po+Pi ■ (5.100) 

Fourier Analyse (drop l's) 

mn(—iu)\ = nqEj — ik.p (5.101) 

The key question: how to relate p to v 
Answer: Equation of state + Continuity 
State 

pn^ 1 = const. =^ (p + pi) (n + ni)~^ = pon 7 (5.102) 
Use Taylor Expansion 



(Po + Pi) ( n o + ni) 7 ~ p n 7 

Hence 



, . Pi ni 

1-| ry 

Po n 



Continuity 



Linearise: 



Fourier Transform 



(5.103) 



^ = 7 ^ (5.104) 
Po n 



— + V. (nv) = (5.105) 
at 

dn-\ On 

— ^ + V. (novi) = =4- — + n V.v = (5.106) 

— ium\ + n ik.V! = (5.107) 
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i.e. 



k.v 



ni = n - 



Combine State & Continuity 



Hi n 

Pi =Pol— = Pol 

n n Q 



k.v 



k.V 



Pol- 



Hence Momentum becomes 



i ■ \ ik Pol, 
mn {— tuj) v = nqh, k.v 



(5.108) 



(5.109) 



(5.110) 



Notice Transverse waves have k.v = 0; so they are unaffected by pressure. 

Therefore we need only consider the longitudinal wave. However, for consistency let us 
proceed as before to get the dielectric tensor etc. 

Choose axes such that k = ke z then obviously: 



v x = E x 

uom 



iq 



E„ 



y uum y 



m —iu + (ik 2, ypo/mnuj) 



Hence 



mq 
com 



1 
1 









fe 2 p 



e = l + 



icr 



e uj 



1-^ 





1 












(5.111) 
(5.112) 

(5.113) 
(5.114) 



1 - ^rpw 

mn 

(Taking account only of 1 species, electrons, for now.) 

We have confirmed the previous comment that the transverse waves (E x , E y ) are unaffected. 
The longitudinal wave is. Notice that e now depends on k as well as oo. This is called 'spatial 
dispersion'. 

For completeness, note that the dielectric tensor can be expressed in general tensor notation 
as 

,2 / 



UJ 2 



1 + kk 



1 - 



1 

k 2 po7 



- 1 



= 1 



1 + kk- 



fc 2 P07 



(5.115) 



This form shows isotropy with respect to the medium: there is no preferred direction in 
space for the wave vector k. 
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But once k is chosen, e is not isotropic. The direction of k becomes a special direction. 
Longitudinal Waves: dispersion relation is 

e zz = (as before) (5.116) 

which is 

oj 2 

1 ^ — = . (5.117) 

mn 



or 

cu 2 = cu 2 + A; 5 

mn 



Recall p — n oT = nT; so this is usually written: 

u 2 = oo 2 p + k 21 — = u 2 p + k 2 -fv 2 (5.119) 



m 



[The appropriate value of 7 to take is 1 dimensional adiabatic i.e. 7 = 3. This seems plausible 
since the electron motion is 1-d (along k) and may be demonstrated more rigorously by kinetic 
theory] 

The above dispersion relation is called the Bohm-Gross formula for electron plasma waves. 
Notice the group velocity: 

duj 1 duo 2 'jkv? 

v 9 = -jt = htt-jt = — L - j — 1 + 0. (5.120) 



dk 2u dk 



(u 2 p + 7/c 2 w 2 ) : 



and for kv t > oj p this tends to 72 u t . In this limit energy travels at the electron thermal 
speed. 

5.4.1 Refractive Index Plot 

Bohm Gross electron plasma waves: 

^ = r4fi-^ (5.121) 



Transverse electromagnetic waves: 



7eV 2 e 1 UJ 2 



N 2 = 1^1 — -JJ (5.122) 

These have just the same shape except the electron plasma waves have much larger vertical 
scale: 

On the E-M wave scale, the plasma wave curve is nearly vertical. In the cold plasma it was 
exactly vertical. 

We have relaxed the Cold Plasma approximation. 
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Figure 5.8: Refractive Index Plot. Top plot on the scale of the Bohm-Gross Plasma waves. 
Bottom plot, on the scale of the E-M transverse waves 



5.4.2 Including the ion response 

As an example of the different things which can occur when ions are allowed to move include 
longitudinal ion response: 

UJ 2 cu 2 

^ tzz ^ U 2 — fc2pe > ^2 _ k 2 pgi (5.123) 
m e n e "'" niiUi 

This is now a quadratic equation for u> 2 so there are two solutions possible for a given oo. One 
will be in the vicinity of the electron plasma wave solution and the inclusion of which is 
<C u 2 e will give a small correction. 

Second solution will be where the third term is same magnitude as second (both will be 
1). This will be at low frequency. So we may write the dispersion relation approximately 

as: 



Ul LU 2 



i.e. 



Lu — 

m e n e mini 



j k 2 Pili ^pi k 2 p, ' < 



= (5.124) 



rriini u 2 e m e n e 
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= k 2 
= k 2 



IJPi JePe 

rii n e 
liTi + 7 e T e 



1 

rrii 



rrii 



(5.125) 



[In this case the electrons have time to stream through the wave in 1 oscillation so they 
tend to be isothermal: i.e. j e — 1. What to take for 7, is less clear, and less important 
because kinetic theory shows that these waves we have just found are strongly damped 
unless Ti < T e ] 

These are 'ion- acoustic 7 or 'ion-sound' waves 



c s is the sound speed 



— c 

k 2 s 



liTi + T e 



rrii 



rrii 



(5.126) 



(5.127) 



Approximately non-dispersive waves with phase velocity c s . 



5.5 Electrostatic Approximation for (Plasma) Waves 

The dispersion relation is written generally as 

N A (N A E) + e.E = N(N.E) - A 2 E + e.E = (5.128) 

Consider E to be expressible as longitudinal and transverse components E £ , E t such that 
N A = 0, N.E t = 0. Then the dispersion relation can be written 

N (N.E £ ) - N 2 (E e + Et) + e. (E* + E t ) = -N 2 E t + e.E t + e.E e = (5.129) 

or 

(A 2 - e) .E t = e.E e (5.130) 

Now the electric field can always be written as the sum of a curl-free component plus a 
divergenceless component, e.g. conventionally 

E= -V0 + v __A_^ (5.131) 

Curl-free Divergence- free 

Electrostatic Electromagnetic 

and these may be termed electrostatic and electromagnetic parts of the field. 

For a plane wave, these two parts are clearly the same as the longitudinal and transverse 
parts because 

— V0 = — ik0 is longitudinal (5.132) 
and if V.A = (because V.A = (w.l.o.g.)) then k.A = so A is transverse. 



115 



'Electrostatic ' waves are those that are describable by the electrostatic part of the electric 
field, which is the longitudinal part: \E e \ 3> \E t \. 

If we simply say E t = then the dispersion relation becomes e.E^ = 0. This is not the most 
general dispersion relation for electrostatic waves. It is too restrictive. In general, there is 
a more significant way in which to get solutions where \Ef \ ^> \E t \. It is for N 2 to be very 
large compared to all the components of e : iV 2 3>|| e ||. 

If this is the case, then the dispersion relation is approximately 

N 2 E t = e.E e ; (5.133) 

E t is small but not zero. 

We can then annihilate the E 4 term by taking the N component of this equation; leaving 

N.e.E* = (N.e.N) E e = : k.e.k = . (5.134) 



When the medium is isotropic there is no relevant difference between the electrostatic dis- 
persion relation: 

N.e.N = (5.135) 

and the purely longitudinal case e.N = 0. If we choose axes such that N is along z, then the 
medium's isotropy ensures the off-diagonal components of e are zero so N.e.N = requires 
tzz = =>- e.N = 0. However if the medium is not isotropic, then even if 

N.e.N (= N 2 e zz ) = (5.136) 

there may be off-diagonal terms of e that make 

e.N ^ (5.137) 

In other words, in an anisotropic medium (for example a magnetized plasma) the electrostatic 
approximation can give waves that have non-zero transverse electric field (of order ||e||/iV 2 
times Eg) even though the waves are describable in terms of a scalar potential. 

To approach this more directly, from Maxwell's equations, applied to a dielectric medium 
of dielectric tensor e, the electrostatic part of the electric field is derived from the electric 
displacement 

V.D = V. (e e.E) = p = (no free charges) (5.138) 
So for plane waves = k.D = k.e.E = ik.e.k0. 

The electric displacement, D, is purely transverse (not zero) but the electric field, E then 
gives rise to an electromagnetic field via V A H = <9D/<9t. If N 2 ^>|| e || then this magnetic 
(inductive) component can be considered as a benign passive coupling to the electrostatic 
wave. 

In summary, the electrostatic dispersion relation is k.e.k = 0, or in coordinates where k is 
in the z-direction, e zz = 0. 
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5.6 Simple Example of MHD Dynamics: Alfven Waves 



Ignore Pressure & Resistance. 



Linearize: 



E + V AB = 

V = Vi, B = B + B x (B uniform), j = jj. 
3V 



at 



j A B c 



Fourier Transform: 



E + V A B = 

p(-iu)V = j A B 
E + V A B = 



Eliminate V by taking 5.144 ABo and substituting from 5.145 

E + — (j A B ) A B = 



or 



E 



-mop 



' -{aBoJBo-SSj}- ^ 



— zcup — zcup 

So conductivity tensor can be written (z in 5 direction). 



-lujp 

R2 
-°0 



1 
1 

oo 



where oo implies that En = (because of Ohm's law). Hence Dielectric Tensor 



e = 1 + 



— iue 



1 + 



1 
1 

oo 



Dispersion tensor in general is: 



D 



L) 



NN — N + e 



(5.139) 
(5.140) 

(5.141) 
(5.142) 
(5.143) 

(5.144) 
(5.145) 

(5.146) 
(5.147) 



(5.148) 



(5.149) 



(5.150) 



Dispersion Relation taking N± = N x , N y — 



IDI 



-^ll+l+eo^ 



N\ Nw 







, N 2_ N 2 ± + 1+ _J>_ 





N±N l{ 


'OO' 



(5.151) 
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Figure 5.9: Compressional Alfven Wave. Works by magnetic pressure (primarily). 



Meaning of oo is that the cofactor must be zero i.e. 



-m + 1 



6 B< 



-N A + 1 + 







(5.152) 



The l's here come from Maxwell displacement current and are usually negligible (N]_ ^> 1). 
So final waves are 



1. N 2 = j-jj2 =>■ Non-dispersive wave with phase and group velocities 



c (c 2 e Q B T 



V * = v *-N 



B^_ 
Pop 



where we call 

= va the 'Alfven Speed' 

.PoP. 

Polarization: 

E ll= E z = 0, E x = 0. E y ^0 ^V y = V x ^0 (V z = 0) 



(5.153) 



(5.154) 



(5.155) 

Party longitudinal (velocity) wave — > Compression "Compressional Alfven Wave". 

k 2 2 

2. N 2 = = %- 

Any has unique k\\. Wave has unique velocity in || direction: va- 
Polarization 



E z = E y = E x ^0 ^V x = Vy^O (V z = 0) 
Transverse velocity: "Shear Alfven Wave". 

Works by field line bending (Tension Force) (no compression). 



(5.156) 
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Figure 5.10: Shear Alfven Wave 



5.7 Non-Uniform Plasmas and wave propagation 

Practical plasmas are not infinite & homogeneous. So how does all this plane wave analysis 
apply practically? 

If the spatial variation of the plasma is slow c.f. the wave length of the wave, then coupling 
to other waves will be small (negligible). 



Figure 5.11: Comparison of sudden and gradualy refractive index change. 

For a given u, slowly varying plasma means N/^ ^> A or kN/^ ^> 1. Locally, the plasma 
appears uniform. 

Even if the coupling is small, so that locally the wave propagates as if in an infinite uniform 
plasma, we still need a way of calculating how the solution propagates from one place to 
the other. This is handled by the 'WKB(J)' or 'eikonal' or 'ray optic' or 'geometric optics' 
approximation. 

WKBJ solution 

Consider the model 1-d wave equation (for field u) 





(5.157) 



with k now a slowly varying function of x. Seek a solution in the form 

E = exp (icf) (x)) {—iujt implied) 



(5.158) 
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is the wave phase (= kx in uniform plasma). 
Differentiate twice 



ax 2 iix 2 \ dx 
Substitute into differential equation to obtain 

/ r i \ 2 



2 , 



S =* (5 ' 160) 



Recognize that in uniform plasma ^ = 0. So in slightly non-uniform, 1st approx is to 
ignore this term. 

^ ~ ±h(x) (5.161) 
ax 

Then obtain a second approximation by substituting 

dx-^^x (5 - 162) 



so 



~ k 2 ±t^ (5.163) 



dx y dx 
ax \ 2A; ax 



~ ± I k ± — — using Taylor expansion. (5.164) 



Integrate: 



Hence E is 



<f) ~ ± J £;dx -M hi (fc 2 ) (5.165) 



E = e i(f> = exp (±i ^ A;dx) (5.166) 

This is classic WKBJ solution. Originally studied by Green & Liouville (1837), the Green 
of Green's functions, the Liouville of Sturm Liouville theory. 

Basic idea of this approach: (1) solve the local dispersion relation as if in infinite homogeneous 
plasma, to get k(x), (2) form approximate solution for all space as above. 

Phase of wave varies as integral of kdx. 

In addition, amplitude varies as This is required to make the total energy flow uniform. 

fc2 



5.8 Two Stream Instability 

An example of waves becoming unstable in a non-equilibrium plasma. Analysis is possible 
using Cold Plasma techniques. 

Consider a plasma with two participating cold species but having different average velocities. 
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These are two "streams" . 

Speciesl Species2 

Moving. Stationary. (5.167) 
Speed v 

We can look at them in different inertial frames, e.g. species (stream) 2 stationary or 1 
stationary (or neither). 

We analyse by obtaining the susceptibility for each species and adding together to get total 
dielectric constant (scalar 1-d if unmagnetized) . 

In a frame of reference in which it is stationary, a stream j has the (Cold Plasma) suscepti- 
bility 

Xi = ^ • (5-168) 
If the stream is moving with velocity Vj (zero order) then its susceptibility is 

— J 1 - 

Xj = — — 2 • (k & v i m same direction) (5.169) 

(u-kvj) 

Proof from equation of motion: 

-^-E = + v.Vv = (-iu + ik.Vj) v = -i(u- kvj) v . (5.170) 

Current density 

j = Pi v i +Pj.v + pvj. (5.171) 

Substitute in 

V.j + ^ = ik.irpj + ik.vp - icup = (5.172) 

Pi = Pj—J777 (5.173) 



LO — k.Vj 



Hence substituting for v in terms of E: 



X ,e V.E = p 3 = m kE (5 . 174) 

TOj [UJ — K.Vj) 



which shows the longitudinal susceptibility is 



2 



Xi = l — = (5-175) 

m,j€o (uu 2 — kvj) (to — kvj) 

Proof by transforming frame of reference: 

Consider Galileean transformation to a frame moving with the stream at velocity v,. 

x = x / + v i t ; t' = t (5.176) 
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exp % (k.x — uot) = exp % (k.x' — (uo — k.v,) i') 
So in frame of the stream, uo' = uo — k.v.,. 
Substitute in stationary cold plasma expression: 



(5.177) 



2 2 

v _ _ c< ^pi ^pj 

Xl " uo' 2 " (uo-kvA 2 ' 



(5.178) 



Thus for n streams we have 



M7 



2 • 



Longitudinal wave dispersion relation is 



e = 0. 



Two streams 



= e = 1 



uo 



uo 



P2 



(uo — kvi) 2 (uo — kv2) 2 



(5.179) 

(5.180) 
(5.181) 



For given real k this is a quartic in uo. It has the form: 




Figure 5.12: Two-stream stability analysis. 

If e crosses zero between the wells, then 3 4 real solutions for uo. (Case B). 
If not, then 2 of the solutions are complex: uo = uo r ± iuoi (Case A). 
The time dependence of these complex roots is 



exp (—iuot) = exp (—iuo r t ± U3{t) . 



(5.182) 
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The +ve sign is growing in time: instability. 

It is straightforward to show that Case A occurs if 



\k(v 2 -vi)\ < + o; p 3 2 



2 2 



3 
■> 



(5.183) 



Small enough k (long enough wavelength) is always unstable. 

Simple interpretation {uj^ 2 <C u^, f i = 0) a tenuous beam in a plasma sees a negative e if 

~ ^pi- 
Negative e implies charge perturbation causes E that enhances itself: charge (spontaneous) 
bunching. 



Wave damping is due to wave-particle resonance. To treat this we need to keep track of the 
particle distribution in velocity space — > kinetic theory. 

5.9.1 Vlasov Equation 

Treat particles as moving in 6-D phase space x position, v velocity. At any instant a particle 
occupies a unique position in phase space (x, v). 

Consider an elemental volume <i 3 x<i 3 v of phase space [dxdydzdv x dv y dv z ], at (x, v). Write 
down an equation that is conservation of particles for this volume 



5.9 Kinetic Theory of Plasma Waves 




= [v x f (x + dxx, v) - v x f (x, v)] dydzd 3 v 

+ same for dy, dz 

+ [a x f (x, v + dv x ±) - a x f (x, v)] d 3 xdv y dv z 

+ same for dv y , dv z 



(5.184) 




Figure 5.13: Difference in flow across x-surfaces (+y + z). 



a is "velocity space motion", i.e. acceleration. 
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Divide through by d 3 x.d 3 v and take limit 

= V.(v/) + V„.(a/) (5.185) 
[Notation: Use £ <-»• V; ^ <— > v., . 

Take this simple continuity equation in phase space and expand: 

^ + (V.v) / + (v.V) / + (V w .a) / + (a.V„) / = 0. (5.186) 
Recognize that V means here ^ etc. keeping v constant so that V.v = by definition. So 

Now we want to couple this equation with Maxwell's equations for the fields, and the Lorentz 
force 

a = — (E + v A B) (5.188) 
m 

Actually we don't want to use the E retaining all the local effects of individual particles. We 
want a smoothed out field. Ensemble averaged E. 

Evaluate 

V„.a = V„.— (E + v A B) = — V„. (v A B) (5.189) 
m m 

= — B. (V„ A v) = 0. (5.190) 

m 

So RHS is zero. However in the use of smoothed out E we have ignored local effect of one 
particle on another due to the graininess. That is collisions. 

Boltzmann Equation: 

9f df df fdf\ 

\ / collisions 

Vlasov Equation = Boltzman Eq without collisions. For electromagnetic forces: 

^ +v .^ + ^( E + vAB )^=0. (5.192) 
ot ox m ov 



Interpretation: 

Distribution function is constant along particle orbit in phase space: J^/ = 0. 



d^ dl dxdl dvdl ,, 1Q ^ 
dt J ~ dt + dfd* dfdv 15 9dj 

Coupled to Vlasov equation for each particle species we have Maxwell's equations. 
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Vlasov-Maxwell Equations 

^ + v.|^ + ^(E + vAB).|^ = (5.194) 
at ax rrij OVj 

VAE = — , VAB = /i j + -— (5.195) 

V.E = - , V.B = (5.196) 
Coupling is completed via charge & current densities. 

P = EW = E%7/A (5-197) 

3 J 

j = Efc n i V i = Efc //j-vrfV- (5-198) 

Describe phenomena in which collisions are not important, keeping track of the (statistically 
averaged) particle distribution function. 

Plasma waves are the most important phenomena covered by the Vlasov-Maxwell equations. 
6-dimensional, nonlinear, time-dependent, integral-differential equations! 

5.9.2 Linearized Wave Solution of Vlasov Equation 

Unmagnetized Plasma 

Linearize the Vlasov Eq by supposing 

/ = / (v) + A(^expi(k.x-^) , /i small. (5.199) 
also E = Eiexpi (k.x - out) B = B 1 exp i (k.x - out) (5.200) 

Zeroth order f equation satisfied by J^, = 0. First order: 

- iuh + v.ik/x + -(E 1 +vAB 1 ).^=0. (5.201) 
m Ov 

[Note v is not per se of any order, it is an independent variable.] 

Solution: 

h = \ ~ (Ei + v A Bi) (5.202) 
i [u — k.v) m ov 

For convenience, assume f is isotropic. Then ^ is in direction v so v A B^^ = 

i(uo — k.v) 

We want to calculate the conductivity cr. Do this by simply integrating: 

j = / qf^dh = — / ^—d z v .Ei. (5.204) 

7 %m J uj — k.v 
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/! = m (5.203) 



Here the electric field has been taken outside the -u-integral but its dot product is with 
dfo/dv. Hence we have the tensor conductivity, 



g2 , „dfo 



[ V a , v d 3 v (5.205) 
J uj — k.v 



vm 

Focus on zz component: 

_ „2 r v dfo 

1 + X ZZ = £ ZZ = 1 + —^ = 1 + — ~^d\ (5.206) 
— tcoeo ujmeo J 00 — k.v 

Such an expression applies for the conductivity (susceptibility) of each species, if more than 
one needs to be considered. 

It looks as if we are there! Just do the integral! 

Now the problem becomes evident. The integrand has a zero in the denominator. At least 
we can do 2 of 3 integrals by defining the 1-dimensional distribution function 

f z (v z ) = J f(v)dv x dv y (k = kz) (5.207) 



r v ^ 

/ Zdv ; dv z (5.208) 
J uj — kv z 



Then 

q 2 r „ 9L 

X = 

(drop the z suffix from now on. 1-d problem). 

How do we integrate through the pole at v — |? Contribution of resonant particles. Crucial 
to get right. 

Path of velocity integration 

First, realize that the solution we have found is not complete. In fact a more general solution 
can be constructed by adding any solution of 

[We are dealing with 1-d Vlasov equation: |£ + + ^ |£ = 0.] Solution of this is 

fi = g(vt-z,v) (5.210) 
where g is an arbitrary function of its arguments. Hence general solution is 

_2_ E— 

fx = — ^ - exp % (kz -ut)+g (vt - z, v) (5.211) 
i{u — kv) 

and g must be determined by initial conditions. In general, if we start up the wave suddenly 
there will be a transient that makes g non-zero. 
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So instead we consider a case of complex to (real k for simplicity) where uo = u r + iuo,i and 
uji > 0. 



This case corresponds to a growing wave: 

exp(— iuit) = exp(—iu r t + Lurf) 



(5.212) 



Then we can take our initial condition to be j\ = at t — > — oo. This is satisfied by taking 
^7 = 0. 

For Ui > i/ie complementary function, g, is zero. 

Physically this can be thought of as treating a case where there is a very gradual, smooth 
start up, so that no transients are generated. 

Thus if Ui > 0, the solution is simply the velocity integral, taken along the real axis, with 
no additional terms. For 



u)i > 0, 



X 



q 



9f 
D Bv 



-dv 



ume Jc oj — kv 

where there is now no difficulty about the integration because u is complex. 









.a 




t k . 















(5.213) 



Figure 5.14: Contour of integration in complex v -plane. 
The pole of the integrand is at v — ¥ which is above the real axis. 

The question then arises as to how to do the calculation if u>i < 0. The answer is by "analytic 
continuation", regarding all quantities as complex. 

"Analytic Continuation" of x is accomplished by allowing u/k to move (e.g. changing the a>j) 
but never allowing any poles to cross the integration contour, as things change continuously. 



Remember (Fig 5.15) 



j> Fdz = ^2 residues x 2ni 



(5.214) 



(Cauchy's theorem) 

Where residues = \im z _> Zk [F(z) / (z—Zk)] at the poles, Zk, of F(z). We can deform the contour 



how we like, provided no poles cross it. Hence contour (Fig 5.16) 
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Figure 5.15: Cauchy's theorem. 



CO: -O 



Figure 5.16: Landau Contour 



We conclude that the integration contour for Ui < is not just along the real v axis. It 
includes the pole also. 

To express our answer in a universal way we use the notation of "Principal Value" of a 
singular integral defined as the average of paths above and below 



P 



f—dv= 1 \f +[ 
J v — vq 2 Uci Jc 2 \ 



-dv 



(5.215) 



4 




Figure 5.17: Two halves of principal value contour. 



Then 



I 2 , 

X = iP 

oome 



v 



-^—dv 
uj — kv 



-2 ' — — 
2 m k 2 dv 



(5.216) 



Second term is half the normal residue term; so it is half of the integral round the pole. 
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Figure 5.18: Contour equivalence. 

Our expression is only short-hand for the (Landau) prescription: 

"Integrate below the pole". (Nautilus). 

Contribution from the pole can be considered to arise from the complementary function 
g(vt — z,v). If g is to be proportional to exp(ikz), then it must be of the form g = exp[ik(z — 
vt)]h(v) where h(v) is an arbitrary function. To get the result previously calculated, the value 
of h(v) must be (for real u) 



h(v) 



TV 



q_l dfo 
mk dv 



8 v- 



(5.217) 



(so that 



h 



iue n 



-vgdv 



Til 



k 2 dv 



urnie r 



■) 



(5.218) 



This Dirac delta function says that the complementary function is limited to particles with 
"exactly" the wave phase speed |. It is the resonant behaviour of these particles and the 
imaginary term they contribute to x that is responsible for wave damping. 

We shall see in a moment, that the standard case will be Ui < 0, so the opposite of the 
prescription u>i > that makes g — 0. Therefore there will generally be a complementary 
function, non-zero, describing resonant effects. We don't have to calculate it explicitly 
because the Landau prescription takes care of it. 



5.9.3 Landau's original approach. (1946) 

Corrected Vlasov's assumption that the correct result was just the principal value of the inte- 
gral. Landau recognized the importance of initial conditions and so used Laplace Transform 
approach to the problem 

_ roc 

A(p) = / e- pt A(t)dt (5.219) 



o 



The Laplace Transform inversion formula is 

1 rs+ioo _ 

A(t) = — / e pt A(p)dp (5.220) 

where the path of integration must be chosen to the right of any poles of A(jp) (i.e. s large 
enough). Such a prescription seems reasonable. If we make dt(p) large enough then the A(p) 
integral will presumably exist. The inversion formula can also be proved rigorously so that 
gives confidence that this is the right approach. 
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If we identify p — > —iuj, then the transform is A — J e tult A{t)dt, which can be identified as 
the Fourier transform that would give component A oc e~ lw *, the wave we are discussing. 
Making dt{p) positive enough to be to the right of all poles is then equivalent to making 5s(uj) 
positive enough so that the path in w-space is above all poles, in particular uji > 5s(kv). For 
real velocity, v, this is precisely the condition uj^ > 0, we adopted before to justify putting 
the complementary function zero. 

Either approach gives the same prescription. It is all bound up with satisfying causality. 



5.9.4 Solution of Dispersion Relation 

We have the dielectric tensor 



e= l+x= 1+ 



ujme 



P 



-dv — ni 



u — kv 



k 2 dv 



(5.221) 



for a general isotropic distribution. We also know that the dispersion relation is 

= (-iV 2 + e t ) 2 e = (5.222) 



-N 2 + e t 
-N 2 + e t 
e 



Giving transverse waves N 2 = e t and 

longitudinal waves e = 0. 

Need to do the integral and hence get e. 

Presumably, if we have done this right, we ought to be able to get back the cold-plasma 
result as an approximation in the appropriate limits, plus some corrections. We previously 
argued that cold-plasma is valid if | ^> v t . So regard ^ as a small quantity and expand: 



P 



,dfo 
dv 



dv 



1 



dv 



1 + 



kv 



kv 

LU 



UJ J 



UJ 



n + 



2kv (kv 

1 + + 3 — 

UJ \uj / 

3nT k 2 

+ ... 



+ .. 



+ 



dv 



dv (by parts) 



m uj' 



(5.223) 



Here we have assumed we are in the particles' average rest frame (no bulk velocity) so that 
/ fovdv = and also we have used the temperature definition 

nT = J mv 2 f dv , (5.224) 

appropriate to one degree of freedom (1-d problem). Ignoring the higher order terms we get: 




T k 2 VI dfo 

3 2 +7 ™7¥~ IT 

muj z k z n ov 



(5.225) 
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This is just what we expected. Cold plasma value was e = 1 — . We have two corrections 

1. To real part of e, correction = 3 due to finite temperature. We could have 
got this from a fluid treatment with pressure. 

2. Imaginary part — > antihermitian part of e — > dissipation. 

Solve the dispersion relation for longitudinal waves e = (again assuming k real u complex). 
Assume Ui <C uj r then 

x2 2 2 r T ^ ^ 2 l<9/o, . 

p mw 2 k 2 n dv k 
9 r T A; 2 w r 2 l%, , 



r 



1 2 1 <9/ 1 2 nu r l df 

HenCeWl " ^Fn^'^S^n^ 1 - (5 ' 227) 



For a Maxwellian distribution 

The difference between cj r and cj p may not be important in the outside but ought to be 
retained inside the exponential since 



m oj 2 
2T¥ 



\ ,Tk 2 
1 + 3 



m uj 2 



moo 2 3 

2T# + 2 (5 ' 23i: 



/7T\ § ^ 1 / mcj 2 3\ 
So ^-^) ^-p(-^--) (5-232) 

Imaginary part of uj is negative =>- damping. This is Landau Damping. 

Note that we have been treating a single species (electrons by implication) but if we need 
more than one we simply add to x- Solution is then more complex. 



5.9.5 Direct Calculation of Collisionless Particle Heating 

(Landau Damping without complex variables!) 

We show by a direct calculation that net energy is transferred to electrons. 
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Suppose there exists a longitudinal wave 

E = E cos(kz - iut)z (5.233) 

Equations of motion of a particle 

^ = —Ecos(kz-out) (5.234) 

ft = V (5 - 235) 

Solve these assuming E is small by a perturbation expansion v — v + Vi + z — z (t) + 
z 1 (t) + ... . 

Zeroth order: 

dv Q 

——- = v o = const , z = Zi + vot (5.236) 

where z^ = const is the initial position. 
First Order 

E cos {kz - out) = —E cos (k (z t + v t) - out) (5.237) 



dvi 


Q_ 


~dt 


m 


dz\ 




~dt 


= Vl 



(5.238) 

at 

Integrate: 

qE sin (kz { + kv - out) 

Vi = : h const. (5.239) 

m kv — oj 

take initial conditions to be i>i,i>2 = 0. Then 

_ qE sin (kzj + Aut) - sin (kzj) (^OAn) 
Vl ~ m Alu ' ^ ' ' 



where Aoj = kv — ou, is (-) the frequency at which the particle feels the wave field. 

qE 



zi = 

m 



cos kzi — cos (kzi + Aut) sin kzi 
Au~ 2 1 Au 



(5.241) 



(using ^(0) = 0). 

2nd Order (Needed to get energy right) 
dv2 qE r . 

— — = ^rjcos [kzi + kv t — out + kz\) — cos [kzi + kv t — out)} 
dt M 

qE 

= — kzi{-sm(kzi + Aut)} {kzx < 1) (5.242) 
Now the gain in kinetic energy of the particle is 

\nv 2 - ^mvl = ^m{( y Vo + v 1 +v 2 + ...f -v%} 

= ^-{2v v 1 + v\ + 2v v 2 + higher order} (5.243) 
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and the rate of increase of K.E. is 



d (\ 2 \ / dv\ dvi dv 2 



dU2 Wj = m ^ _ + „,_+„„_ j (5.244) 

We need to average this over space, i.e. over Zi. This will cancel any component that simply 
oscillates with Zj. 

d fl 2 \\ I dv\ dv\ dv 2 y 



dt \2 mv ) =\ Vo ^ +Vi ^ +Vo ^ r (5 - 245) 



Vl 



dvi 
dt 
dvi 
dt 



v - 



dv 2 
~dt 



(5.246) 



'q 2 E 2 
m 2 



sin {kzi + Aut) — sin kzi 
Au 



cos (kzi + Aut) 



q 2 E 2 I sin (kzi + Aut) — sin (kzi + Aut) cos Aut + cos (kzi + Aut) sin Aut 



to" 



Au 



cos (kzi + Aut) ^ 

o 2 E 2 /sin Aut 9 ,, . . \ 
cos 2 (fczj + Acut) \ 



to" \ Ao; 
g 2 £ 2 1 sin Aut 
m 2 2 Ac<j 



-q 2 E 2 



-kvi 



o 



cos kz;, — cos (fc^i + Aut) sin \ . , , , y 

- t— — - I sin (kzi + Aut) 



-q 2 E 2 

m 2 
q 2 E 2 kv 



kvn 



\ Au 2 

/ sin Aut cos Aut \ 
-t 



to" 



V Au 2 
sin Aut 



Au 2 



+ t- 



Au J 
cos Aut 



Au 



sin 2 (kzi + Aut) 



Au 



(5.247) 



(5.248) 



Hence 



d 1 

dt2 



mv 



q 2 E 2 
2m 

q 2 E 2 r 



sin Aa>t 
Au; 



sin Acut 



2m 



-u sin Ac<jt 



Aw 2 



+ 



Au 2 
ut 



kv t 



cos Acut 



Au 



Au 

cos Acjt + t cos Acut 



(5.249) 
(5.250) 



This is the space-averaged power into particles of a specific velocity vq. We need to integrate 
over the distribution function. A trick identify helps: 



—u . ut 

— — - sin Aut + — — cos Aut + t cos Aut 
Au 2 Au 



d ( u sin Aut 



dAu 



Au 



1 J}_ 
k Ovq 



u sin Aut 
Au 



+ sin Aut 



+ sinAcutj (5.251) 
(5.252) 
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Hence power per unit volume is 



d 1 



dt2 
q 2 E 2 



m v 2 ) f (vq) dv 



2mk 
q 2 E 2 
2mk 



f(v 



d (ou sin Auut 



dvo \ Acu 
ou sin Aout 



+ sin Aout J dv 
df 



Aou 



+ sin Aout ) —^dvo 
J dv 



(5.253) 



As t becomes large, sinAwt = sin(/ct> — cj)t becomes a rapidly oscillating function of v . 
Hence second term of integrand contributes negligibly and the first term, 



oj sin Aut sin Aut 



oc 



Au 



Aut 



cut 



(5.254) 



becomes a highly localized, delta-function-like quantity. That enables the rest of the inte- 
grand to be evaluated just where Acu = (i.e. kvo — cu = 0). 




Figure 5.19: Localized integrand function. 



So: 



x = Acut 
and / s -^dz 

J X 



(kv — cu)t. 

71 SO 



q 2 E 2 cu df 
2mk k dv 



I 



SIM 

X 



dx 



P = -E 



7rq 2 cu df . 



(5.255) 



(5.256) 



2mk 2 dv ¥ 

We have shown that there is a net transfer of energy to particles at the resonant velocity 
from the wave. (Positive if |£| is negative.) 



5.9.6 Physical Picture 

Acu is the frequency in the particles' (unperturbed) frame of reference, or equivalently it is 
kv' where v' is particle speed in wave frame of reference. The latter is easier to deal with. 
Acut = kv' t is the phase the particle travels in time t. We found that the energy gain was of 
the form 

sin Acut , , . . . 

-d(Acut). (5.257) 



/ 



Acut 
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I-* Oo-t *1 




Figure 5.20: Phase distance traveled in time t. 

This integrand becomes small (and oscillatory) for Aut ^> 1. Physically, this means that 
if particle moves through many wavelengths its energy gain is small. Dominant contribution 
is from Aut < it. These are particles that move through less than \ wavelength during the 
period under consideration. These are the resonant particles. 




Figure 5.21: Dominant contribution 
Particles moving slightly faster than wave are slowed down. This is a second-order effect. 




Figure 5.22: Particles moving slightly faster than the wave. 

Some particles of this vq group are being accelerated (A) some slowed (B). Because A's are 
then going faster, they spend less time in the 'down' region. B's are slowed; they spend more 
time in up region. Net effect: tendency for particle to move its speed toward that of wave. 

Particles moving slightly slower than wave are speeded up. (Same argument). But this is 
only true for particles that have "caught the wave" . 

Summary: Resonant particles' velocity is drawn toward the wave phase velocity. 

Is there net energy when we average both slower and faster particles? Depends which type 
has most. 

Our Complex variables wave treatment and our direct particle energy calculation give con- 
sistent answers. To show this we need to show energy conservation. Energy density of 
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Figure 5.23: Damping or growth depends on distribution slope 



wave: 



W 



1 1 |r o, I 
_2^ 2_!!_, , y 

<sin 2 > Electrostatic Particle Kinetic 



+ n-m\v 2 \ 
2 l l 



(5.258) 



Magnetic wave energy zero (negligible) for a longitudinal wave. We showed in Cold Plasma 
treatment that the velocity due to the wave is v — qE Hence 

,2" 

(again electrons only) (5.259) 



W 



l e E 2 
2 2 



UJ 2 



(5.260) 



When the wave is damped, it has imaginary part of uj, u>i and 

dW — 1 dE 2 

= W 

dt E 2 dt 

Conservation of energy requires that this equal minus the particle energy gain rate, P. Hence 

2 7T UJ 1 dfo 2 



_p — I — e 2 nq \ ^ 1 » 

2mk 2 dv I fc 

_ _ 



eoEl 
1 



1 + 



UJ, 



p 2 k 2 n dv 



x 



1 + 



(5.261) 



So for waves such that uj ~ uj p , which is the dispersion relation to lowest order, we get 

2 nuj r 1 dfo 



UJj = UJ. 



p 2k 2 n dv 



(5.262) 



This exactly agrees with the damping calculated from the complex dispersion relation using 
the Vlasov equation. 

This is the Landau damping calculation for longitudinal waves in a (magnetic) field-free 
plasma. Strictly, just for electron plasma waves. 

How does this apply to the general magnetized plasma case with multiple species? 

Doing a complete evaluation of the dielectric tensor using kinetic theory is feasible but very 
heavy algebra. Our direct intuitive calculation gives the correct answer more directly. 
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5.9.7 Damping Mechanisms 

Cold plasma dielectric tensor is Hermitian. [Complex conjugate*, transpose T = original 
matrix.] This means no damping (dissipation). 

The proof of this fact is simple but instructive. Rate of doing work on plasma per unit 
volume is P = E.j. However we need to observe notation. 

Notation is that E(k, u) is amplitude of wave which is really 9?(E(k, u) expi(k.x — ut)) and 
similarly for j. Whenever products are taken: must take real part first. So 



3? (Eexpi (k.x - out)) M (jexpi (k.x - ut)) 



Ee^ + EV 



-i<f> 



k.x — cut.) 



E.je 2 ^ + E.j* + E*.j + E*.j*e" 2 ^ 



(5.263) 



The terms e 2 ^ & e~ 2i<?i are rapidly varying. We usually average over at least a period. These 
average to zero. Hence 

(P> = i[Ej* + E*j] = i»(Ej*) (5.264) 



Now recognize that j = cr.E and substitute 

1 



4 



E.<r*.E* + E*.<r.El 



But for arbitrary matrices and vectors: 

A. MB = B.M T .A; 
(in our dyadic notation we don't explicitly indicate transposes of vectors). So 



E.<r*.E* = E*.<7* T .E 
1 



hence 



(P) 



-E*. 



<T + (7 



.E 



(5.265) 



(5.266) 



(5.267) 



(5.268) 



If e = 1 H r- — a is hermitian e* T = e, then the conductivity tensor is antihermitian 



<t (if uj is real). In that case, equation 5.268 shows that < P >= 0. No dissipation. 



Any dissipation of wave energy is associated with an antihermitian part of cr and hence e. 
Cold Plasma has none. 



Collisions introduce damping. Can be included in equation of motion 

m— = q (E + v A B) — mv v 
at 



(5.269) 



where v is the collision frequency. 



137 



Whole calculation can be followed through replacing m(—iu) with m{y — iuj) everywhere. 
This introduces complex quantity in S, D, P. 

We shall not bother with this because in fusion plasmas collisional damping is usually neg- 
ligible. See this physically by saying that transit time of a wave is 



Size 



1 meter 



~ 3 x 10 seconds. 



Speed 3 x 10+ 8 m/s 
(Collision frequency) -1 ~ lOfis — > 1ms, depending on T e ,n 



(5.270) 



When is the conductivity tensor Antihermitian? 



Cold Plasma: 



S 
iD 




-iD 
S 







p 



where 



S = 
D 



E 



ELU . 
P3 

3 ~i L0 2 -Q, 2 
, ,2 



pi 

3 u) 2 



(5.271) 



This is manifestly Hermitian if uj is real, and then a is anti-Hermitian. 

This observation is sufficient to show that if the plasma is driven with a steady wave, there 
is no damping, and k does not acquire a complex part. 

Two stream Instability 



= i-E 



P3 



(5.272) 



In this case, the relevant component is Hermitian (i.e. real) if both uj and k are real. 

But that just begs the question: If uj and k are real, then there's no damping by definition. 

So we can't necessarily detect damping or growth just by inspecting the dieletric tensor form 
when it depends on both uj and k. 

Electrostatic Waves in general have e = which is Hermitian. So really it is not enough to 
deal with e or x- We need to deal with a = —iuje x, which indeed has a Hermitian component 
for the two-stream instability (even though x is Hermitian) because uj is complex. 



5.9.8 Ion Acoustic Waves and Landau Damping 

We previously derived ion acoustic waves based on fluid treatment giving 



= 1 



UJ. 



pe 



UJ 



pi 



UJ Z 



m e n e 



UJ* 



mini 



(5.273) 



Leading to uj 2 ~ k 2 



rrii 
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Kinetic treatment adds the extra ingredient of Landau Damping. Vlasov plasma, unmagne- 
tized: 



k 2 Jc v — f dv n 

k 



k 2 



dfd dv 



c v — | dv n 



(5.274) 



Both electron and ion damping need to be considered as possibly important. 

Based on our fluid treatment we know these waves will have small phase velocity relative to 
electron thermal speed. Also c s is somewhat larger than the ion thermal speed. 

















\u . 







Figure 5.24: Distribution functions of ions and electrons near the sound wave speed. 
So we adopt approximations 



Vte > p V tl < (<)- 



(5.275) 



and expand in opposite ways. 
Ions are in the standard limit, so 



V ■ ~ — 

At — 9 



m lo 



\w/k 



(5.276) 



Electrons: we regard f as small and write 



P 



1 df oe dv 
v — | dv n 



- P 



1 df oe dv 



v dv n 
2 rdf 

dv 



- -I 

n J 



n J dv 2 
2 f m 



2T f 



f oe dv for Maxwellian. 



me 



(5.277) 



Write F = f /n. 

Contribution from the pole is as usual so 

,2 



LO. 



Xe k2 



m e . dF 0( 



ui/k 



(5.278) 
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Collecting real and imaginary parts (at real a;) 



e r {uj r ) = 1 + 



UJpeTrie w. 

k 2 T P 



pi 



1 + 



3T, k 2 



m oof. 



— 71 



CO 



2 9F e 



pe 



dv 



u/k + W 



pi 



dv 



\u)/k 



(5.279) 
(5.280) 



The real part is essentially the same as before. The extra Bohm Gross term in ions appeared 
previously in the denominator as 



i .2 , ,2 

2 _ k 2 pj"ti ,,,2 



1 + 



3T,- k 2 



rrij uj a 



(5.281) 



Since our kinetic form is based on a rather inaccurate Taylor expansion, it is not clear that 
it is a better approx. We are probably better off using 



oj 2 1 - 



1 

3Tjk 2 ' 
rriiLO 2 



(5.282) 



Then the solution of e r (uj r ) = is 



k 2 



3T,; 



rrii 



1 + k 2 \ 2 



(5.283) 



De 



as before, but we've proved that 7 e = 1 is the correct choice, and kept the k 2 \ 2 De term (1st 
term of e r ). 

The imaginary part of e gives damping. 



General way to solve for damping when small 

We want to solve e(k, uj) — with u = u r + iu>i , oji small. 
Taylor expand e about real oj t : 

de 

s(uj) ~ e(uj r ) + ivi — \ Ur 

= e(uj r ) + ioJiJ^—e(uj r ) 

OUJ r 

Let uj r be the solution of e r {uj r ) = 0; then 



d 

s{uj) = iEi(uu r ) + iuji- — e(u r ). 

UUJ r 



This is equal to zero when 



de(u> r ) 

duir 



(5.284) 
(5.285) 

(5.286) 
(5.287) 
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If, by presumption, <C e r , or more precisely (in the vicinity of e — 0), dei/duj r <C 
then this can be written to lowest order: 



0>; 



g» (^r) 
du) r 



Apply to ion acoustic waves: 



de r {u r ) 
dco r 



CO, 
CO, 



2 r 

pi 
3 

r L 



2 + 4 



m, ; CO, 



so 



7T Ci£ 

k 2 co 2 pi 



2 + 4^4 



a;, 



pe 



dv 



■j/k + v 



dF n . 



pi 



dv 



\u)/k 



For Maxwellian distributions, using our previous value for u r , 

( m e V 
" V27rT e y 



dv k 



2 m e f _msil 

e 2T <= 



T e + 3T, 



m P 1 + 



3Ti 



cxp 



^1 + PA 2 , V 2m i 1 + k2X D i 





1 H 


^ _ 


2 


T e 


/l- 


\- k2X De 



V27T Vmj 

where the exponent is of order m e /mi here, and so the exponential is 1. And 



dF m i 
dv 



1 rrii 



1 ~k 



1 + 



3Ti 



cxp 



T fi 1 



3T; 



2T, 1 + k 2 X 2 D 



Hence 



COj, 

uo r 




T e 1 + 



37; 



2Ti 1 + PA|, 



T P 1 



1 3 

m e \ 2 

mj \Ti 



cxp 



32; 



2T t 1 + PA 2 - 



electron 



ion damping 



de r /dui r 
(5.288) 



(5.289) 



(5.290) 



(5.291) 



(5.292) 



(5.293) 



(5.294) 
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[Note: the coefficient on the first line of equation 5.294 for u>i/u r reduces to 
Ti/T e < 1 and k\ D e < 1.] 

Electron Landau damping of ion acoustic waves is rather small: ^ r 
Ion Landau damping is large, ~ 1 unless the term in the exponent is large. That is 



vr/8 for 



/ m e 1 

70' 



unless 



Ti 



> 1 



(5.295) 



Physics is that large |? pulls the phase velocity of the wave: y e ^ Ti = c s above the ion 
thermal velocity v t i = \J~^- If c s ^> v t % there are few resonant ions to damp the wave. 

[Note. Many texts drop terms of order |f early in the treatment, but that is not really 
accurate. We have kept the first order, giving extra coefficient 



1 + 



3T,, 



3T\ 
2T P 



(5.296) 



and an extra factor 1 + ^ in the exponent. When Tj ~ T e we ought really to use full 
solutions based on the Plasma Dispersion Function.] 

5.9.9 Alternative expressions of Dielectric Tensor Elements 

This subsection gives some useful algebraic relationships that enable one to transform to 
different expressions sometimes encountered. 



Xz 



q 2 u 



uome Jc u> — kv uo z me k Jc \oo — kv J ov 



q 



2 l 



1 df 
me k 2 Jc t — v dv 



dv 



1 ldf 



k 2 Jc t — v ti Ov 

k 



dv 



CO 



P 

k 2 



P 



1 dF ^ .dF . 

dv — ixi— — W 

— v 2v ov k 



(5.297) 
(5.298) 

(5.299) 

(5.300) 



where F = ^ is the normalized distribution function. Other elements of x involve integrals 
of the form 



come c 



3 dv t 



d 3 v . 



(5.301) 



q* J co — k.v 

When k is in z-direction, k.v = k z v z . (Multi dimensional distribution / ). 

If (e.g., Xxy) I 7^ z an d j 7^ I then the integral over v\ yields J ^dvi = 0. If j — I ^ z then 

dfo 



Jv^dv^-Jfjv, 



(5.302) 



142 



by parts. So, recalling the definition f z = J fdv x dv y , 

Jo 



Xxx Xyy 



Q 



ujme J uj — k.v 



-dv. 



< [ 

uj J 



OZ j 

- — dv z . 
uj J uj — k.v 



(5.303) 



The fourth type of element is 



X.i 



ujme, 



r v ^ 
f Ux dv z 

J uj — k,v 



d 3 v . 



(5.304) 



This is not zero unless f Q is isotropic (= f (v))- 
If / is isotropic 



dfo = dfodv_ = v z df 
dv z dv dv z v dv 



(5.305) 



Then 



r v x 
J uj — 



dfo 

^d 3 v 



d s v 



f VxVz 1 

J uj — k z v z v 

/ 1 7T~ d v = G 

J uj — k z v z ov x 



(5.306) 



(since the ^-integral of df /dv x is zero). Hence for isotropic F D = f /n, with k in the 
z-direction, 



X 



Foz 



dv. 



C ui—kv z z 








(,) JC i.) — hi)-. % 



0+ 
0+ 



^ k ui—kv z dv z ^Vz 



(5.307) 



(and the terms + are the ones that need isotropy to make them zero). 



e t 
e t 
e l 



(5.308) 



where 



uj: 



uj Jc uj — kv 



-dv. 



1 dF 



k 2 Jc v — r dv. 



dv, 



(5.309) 
(5.310) 



k ^ Uz 

All integrals are along the Landau contour, passing below the pole. 
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5.9.10 Electromagnetic Waves in unmagnetized Vlasov Plasma 

For transverse waves the dispersion relation is 

k 2 c 2 „ r9 , f fozdv z 

— = N 2 = e t = 1 - / 5.311 
c<j z uj n Jc {uj — k z v z ) 

This has, in principle, a contribution from the pole at uj — kv z = 0. However, for a non- 
relativistic plasma, thermal velocity is <C c and the EM wave has phase velocity ~ c. Con- 
sequently, for all velocities v z for which f oz is non-zero kv z <^ u. We have seen with the cold 
plasma treatment that the wave phase velocity is actually greater than c. Therefore a proper 
relativistic distribution function will have no particles at all in resonance with the wave. 

Therefore: 



1. The imaginary part of e t from the pole is negligible. And relativisitically zero. 
2. 



€t — 1 

= 1 

~ 1 

~ 1 



-i-r f-( 



<4 1 

CO 2 
UJ 2 



kv z k 2 v 2 
1 + + — ^ + .. 



dv. 



UJ 



UJ C 



k 2 T 



2 — + - 

uj z m 



1 + 



k 2 v 2 

UJ 2 



(5.312) 



Thermal correction to the refractive index iV is small because — v 4- <^ 1. 



Electromagnetic waves are hardly affected by Kinetic Theory treatment in unmagnetized 
plasma. Cold Plasma treatment is generally good enough. 
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